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ABSTRACT 


Social  network  analysis  frequently  uses  the  idea  of  a  clique  in  a  network  to 
identify  key  subgroups  of  highly-connected  members  of  the  network.  We  formulate  the 
maximum  clique  problem  on  undirected  graphs  and  develop  two  algorithms  to  solve  it:  a 
pruning  algorithm  and  an  enumeration  algorithm.  The  pruning  algorithm  successively 
improves  an  upper  bound  on  the  clique  number  of  a  graph,  and  the  enumeration 
algorithm  successively  finds  larger  and  larger  cliques  in  the  graph.  Both  terminate  with  a 
maximum  clique  in  the  graph,  and,  when  run  together,  provide  an  interval  of  uncertainty 
on  the  size  of  a  maximum  clique  in  a  graph  that  converges  to  zero.  We  apply  our 
algorithms  to  real  examples  in  the  modeling  of  terrorist  social  networks,  and  detennine 
that  our  algorithms  are  efficient  and  practical  for  problems  of  moderate  size. 
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EXECUTIVE  SUMMARY 


An  undirected,  simple  graph,  G=(  V,E),  is  defined  by  a  set  of  vertices,  V,  and  a  set 
of  edges,  E,  consisting  of  unordered  pairs  of  distinct  vertices  that  represent  symmetric, 
pairwise  relationships,  or  adjacencies,  between  those  vertices.  A  clique  in  a  graph  G  is  a 
subset  of  vertices,  Ccf,  such  that  every  pair  of  vertices  in  C  corresponds  to  an  edge  in 
E.  Namely,  every  pair  of  vertices  in  C  must  be  adjacent  in  G,  and  so  C  induces  a 
complete,  or  maximally  connected,  subgraph  of  G.  Finding  cliques  of  maximum 
cardinality  in  a  graph  is  a  long-standing  problem  in  graph  theory,  and  is  referred  to  as  the 
maximum  clique  problem,  or  MCP. 

The  MCP  has  important  applications  in  many  different  domains.  Some  examples 
of  military  applications  are  in  Cryptography  and  Cryptanalysis,  Telecommunications, 
particularly  in  Wireless  Networks  and  in  Radio  Frequency  Assignment,  and  Social 
Networks  Analysis.  Since  the  clique  is  considered  as  the  foundational  idea  for  studying 
cohesive  subgroups  in  social  networks,  this  thesis  considers  the  maximum  clique  problem 
with  an  emphasis  on  its  application  to  military  problems,  especially  in  Social  Networks 
Analysis. 

The  primary  question  that  this  thesis  addresses  is  the  following: 

Can  we  identify  and  implement  an  exact  algorithm  to  solve  the  Maximum  Clique 

Problem  (MCP)  on  undirected  graphs,  in  a  reasonable  time  frame? 

To  achieve  the  above  objective  we  selected  two  algorithms  from  the  current 
literature  and  modified  them  to  improve  their  perfonnance,  developing  the  pruning  and 
enumeration  algorithms.  We  then  implemented  both  of  these  in  a  modem,  powerful  and 
widely  used  programming  language,  Java.  Our  testing  involved  applying  both  algorithms 
to  real-world  situations  which  could  be  modeled  by  undirected  graphs:  terrorist  social 
networks. 

We  verified  that  both  algorithms  solve  the  MCP  on  undirected  graphs,  and  are 
quick  on  relatively  small  graphs.  Furthermore,  the  pruning  algorithm  immediately 


establishes  an  upper  bound  on  the  clique  number  of  a  graph  and  then  successively 
improves  this  bound;  thus  the  pruning  algorithm  can  be  tenninated  early  with  a  valid 
upper  bound  on  the  clique  number.  Similarly,  the  enumeration  algorithm  quickly  finds 
small  cliques,  and  successively  discovers  larger  and  larger  cliques  in  the  graph  as  it 
progresses,  each  of  which  provides  a  lower  bound  on  the  clique  number  of  the  graph.  It 
can  be  terminated  early  with  a  valid  lower  bound  on  the  clique  number. 

If  both  are  run  simultaneously  each  provides  a  bound  the  other  cannot,  and  an 
interval  of  uncertainty  can  be  established  that  will  eventually  be  reduced  to  zero,  at  which 
point  a  maximum  clique  will  have  been  found. 

Moreover  analyzing  a  social  network  using  the  Social  Networks  Analysis 
methods  and  measures  maximum  cliques  and  vertices  degree  may  provide  enough 
information  about  the  social  structure  of  the  network  under  investigation.  More 
particularly,  the  maximum  clique(s)  and  the  clique(s)  as  well,  may  give  us  the  most 
cohesive  subgroups  while  the  vertices  with  the  largest  vertex-degrees  may  show  the  most 
central  actors  in  a  social  network.  However,  since  cliques  have  been  criticized  for  their 
restrictive  nature  we  may  consider  and  study  one  of  the  branches  of  the  relaxation  clique 
problem  and  analyze  a  social  network  under  this  concept.  This  may  be  the  subject  for  an 
extension  of  this  current  thesis  and  an  area  for  further  investigation.  It  is  worthwhile  to 
point  out  that  each  algorithm  can  be  modified  in  fairly  straightforward  ways  to  allow 
various  relaxations  of  the  definition  of  a  clique. 
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I.  INTRODUCTION 


It  has  been  said  that  figures  rule  the  world.  Maybe.  But  I  am  sure  that 
figures  show  us  whether  it  is  being  ruled  well  or  badly. 

Goethe 


A.  CLIQUES  IN  SOCIAL  NETWORK  ANALYSIS 

An  undirected,  simple  graph,  G=(V,E),  is  defined  by  a  set  of  vertices,  V,  and  a  set 
of  edges,  E,  consisting  of  unordered  pairs  of  distinct  vertices  that  represent  symmetric, 
pairwise  relationships,  or  adjacencies,  between  those  vertices.  A  clique  in  a  graph  G  is  a 
subset  of  vertices,  Ccf,  such  that  every  pair  of  vertices  in  C  corresponds  to  an  edge  in 
E.  Namely,  every  pair  of  vertices  in  C  must  be  adjacent  in  G,  and  so  C  induces  a 
complete,  or  maximally  connected,  subgraph  of  G.  Finding  cliques  of  maximum 
cardinality  in  a  graph  is  a  long-standing  problem  in  graph  theory,  and  is  referred  to  as  the 
maximum  clique  problem,  or  MCP. 

The  MCP  has  important  applications  in  many  different  domains.  Some  examples 
of  military  applications  are  in  Cryptography  and  Cryptanalysis,  Telecommunications, 
particularly  in  Wireless  Networks  and  in  Radio  Frequency  Assignment,  and  Social 
Networks  Analysis.  In  fact,  according  to  Wassennan  and  Faust  [10],  “The  clique  is  the 
foundational  idea  for  studying  cohesive  subgroups  in  social  networks....”  This  thesis 
considers  the  maximum  clique  problem  with  an  emphasis  on  its  application  to  military 
problems,  especially  in  Social  Networks  Analysis. 

B.  THESIS  OBJECTIVE  AND  RESEARCH  QUESTION 

The  primary  question  that  this  thesis  addresses  is  the  following: 

Can  we  identify  and  implement  an  exact  algorithm  to  solve  the  Maximum  Clique 

Problem  (MCP)  on  undirected  graphs,  in  a  reasonable  time  frame? 

To  achieve  the  above  objective  we  first  review  the  current  literature  on  MCP,  and 
then  select  two  existing  algorithms  related  to  this  problem.  We  implement  efficient 
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versions  of  both  and  apply  them  to  real-world  situations.  We  focus  on  social  networks, 
where  we  find  a  wide  variety  of  important  applications.  We  will  also  consider  interesting 
variations  to  the  initial  problem,  such  as  various  relaxations  of  the  clique  requirement  that 
overcome  objections  to  the  MCP’s  restrictive  nature  and  modeling  disadvantages. 
Finally,  we  analyze  the  behavior  of  our  two  algorithms  and  draw  conclusions  on  their 
performance,  utility  and  effectiveness.  Our  main  intention  is  to  improve  their 
performance  and  obtain  a  new  and  more  efficient  algorithms  that  can  be  applied  to 
military  applications  of  the  maximum  clique  problem  (and  its  relaxations)  for  problems 
of  reasonable  size. 

C.  LITERATURE  REVIEW  AND  ITS  CONNECTION  TO  CURRENT  WORK 

The  decision  version  of  the  MCP,  in  which  the  goal  is  to  decide  whether  a  given 
graph  contains  a  clique  of  a  given  cardinality,  is  one  of  the  first  problems  that  have  been 
proved  to  be  NP-complete  (see  [1]  for  history  and  discussion),  and  therefore  the  MCP  is 
an  NP-Hard  problem. 

Tarjan  [7]  was  one  of  the  first  who  addressed  the  MCP  and  presented  one  of  the 
first  algorithms  that  gave  a  reasonable  and  effective  solution  to  the  problem  at  hand. 
Tarjan  provided  an  algorithm  with  a  running  time  bounded  by  0{n2"  j ,  where  n  is  the 

order  of  a  given  graph  G .  The  basic  algorithm  examines  every  subset  of  the  vertex  set 
V  (  G  )  of  a  given  a  graph  G  with  n  vertices.  The  algorithm  determines  all  subsets  that 
are  cliques,  and  chooses  the  largest  clique  found  as  the  maximum  clique.  Since  the 
number  of  subsets  of  any  set  with  n  elements  is  2" ,  and  it  takes  0{n )  time  to  check  if  a 
subset  fonns  a  clique,  it  follows  that  the  time  upper  bound  for  this  simple  algorithm  to 
solve  the  MCP  is(9(«2") . 

Tarjan  later  succeeded  in  improving  this  basic  algorithm  and  discovered  other, 
improved  algorithms.  The  first  one  had  a  worst-case  time  bound  of  k  (1.286)”  for  some 
constant  k  .  Hence,  within  a  fixed  amount  of  time  this  improved  algorithm  could  analyze 
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a  larger  graph  than  the  basic  algorithm.  A  few  years  later,  Tarjan  and  Trojanowski  [13] 
presented  a  recursive  algorithm  which  detennines  a  maximum  independent  set  of  n  - 

vertex  graph  in  0(2n ^  time  (see  also  [12]). 

Since  the  early  1970s,  many  papers  have  been  published  with  algorithms  for  the 
MCP  [2].  According  to  L.  Babel  in  [3],  earlier  work  diverged  into  two  directions.  The 
first  concerned  algorithms  solving  the  problem  for  arbitrary  graphs  in  exponential  time, 
the  other  restricted  to  special  classes  of  graphs  where  polynomial  methods  could  be 
found. 

Between  September  1992  and  September  1993,  the  Second  DIMACS 
Implementation  Challenge  took  place.  The  purpose  of  the  challenge  was  to  encourage 
high-quality  research  on  empirical  issues  in  combinatorial  optimization.  The  problem  of 
finding  cliques  in  graphs  was  one  of  the  three  problem  classes  that  was  discussed. 
According  to  the  analytical  results  of  this  challenge,  which  were  presented  in  [4],  it  seems 
unlikely  to  have  a  fast,  i.e.,  polynomial  time,  algorithm  to  solve  this  kind  of  problem 
exactly.  Even  finding  an  approximate  solution  quickly  is  improbable  due  to  the  fact  that  it 
is  an  NP-complete  problem.  Moreover,  all  the  papers  presented  during  this  challenge  on 
finding  cliques  in  a  graph  were  a  mixture  of  exact  and  heuristic  methods,  that  is  a  mixture 
of  exact  and  approximate  methods  [4].  We  present  a  brief  description  of  these  two 
algorithm  categories. 

The  exact  algorithms  are  those  that  have  been  proven  mathematically  to  provide 
an  optimal  solution.  Branch-and-bound,  for  example,  is  a  finite  computing  time  method 
that  has  been  widely  implemented  in  the  efforts  of  solving  the  MCP  or  a  part  of  it  [2,  4]. 
Besides  the  quality  of  the  solution  that  this  method  guarantees,  it  also  provides  the 
solution  in  an  acceptable  finite  computing  time. 

On  the  other  hand,  heuristics  algorithms  are  those  that  cannot  guarantee  any 
solution  quality.  Greedy  algorithms,  Neighborhood  search,  and  Tabu  search  are  typical 
examples  of  heuristic  methods  which  also  have  been  considered  in  approaching  the  MCP 
[4,  9], 
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Early  algorithms  for  solving  the  MCP  as  well  as  recent  approaches  included  the 
branch-and-bound  method.  In  [12],  Wood  referred  to  many  fundamental  approaches  to 
the  maximum  clique  problem  which  include  branch-and-bound  algorithms.  First  of  all, 
Wood  considered,  in  this  paper,  the  significance  of  determining  an  upper  and  a  lower 
bound  for  problems  that  are  NP-complete  like  the  MCP.  He  also  presented  a  branch-and- 
bound  algorithm  for  finding  a  maximum  clique  in  a  graph,  distinguishing  two  algorithms 
which  were  the  most  efficient  ones  known  for  the  maximum  clique  problem  until  1997. 
One  of  these  algorithms  was  developed  by  Babel  and  presented  in  [3]  while  the  other 
belongs  to  Balas  and  Xue  and  was  presented  in  [14].  These  algorithms  calculated  lower 
and  upper  bounds,  which  seems  to  be  a  great  method  for  NP-complete  problems. 
Moreover,  according  to  Wood  [12],  Pardalos  and  Xue  identified  in  their  survey  paper 
[15]  the  following  three  key  questions  that  have  arisen  in  a  branch-and-bound  algorithm 
for  the  MCP.  We  quote  from  [12]: 

1 .  How  to  find  a  good  lower  bound,  i.e.,  a  clique  of  large  size? 

2.  How  to  find  a  good  upper  bound  on  the  size  of  maximum  clique? 

3.  How  to  branch,  i.e.,  break  a  problem  into  smaller  subproblems? 

And  so,  due  to  the  hard  nature  of  finding  even  an  approximate  solution  to  the 
MCP  quickly,  we  selected  an  approach  to  the  problem  that  combines  two  separate 
algorithms,  a  pruning  algorithm,  and  an  enumeration  algorithms.  The  first  finds  a 
sequence  of  improved  upper  bounds  on  the  maximum  clique  size,  until  it  finds  a  clique  of 
size  equal  to  the  current  bound,  and  the  latter  uses  backtracking  search  to  build  a 
sequence  of  larger  and  larger  cliques,  and  therefore  an  increasing  sequence  of  lower 
bounds  on  the  maximum  clique  size;  both  find  all  maximum-size  cliques  in  any  given 
graph. 
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II.  FORMULATING  THE  MAXIMUM  CLIQUE  PROBLEM 


Thanks  to  Euler,  Graph  Theory  is  thriving. 
Year  by  year  it  flourishes  and  blossoms, 
Fertilizing  much  of  mathematics 
And  so  rich  in  all  its  applications. 


Bohdan  Zelinka, 
The  graph  theory  hymns 


A.  BOUNDS  AND  ASSUMPTIONS 


In  order  to  develop  upper  and  lower  bounds  on  clique  sizes  in  a  given  network, 
and  to  justify  our  algorithms  in  Chapter  III,  we  establish  here  a  few  key  observations 
concerning  cliques  in  simple,  undirected  graphs. 


1.  Maximum  Clique  Size  Bounds 

In  [8],  it  was  shown  that  every  graph  G  with  n  vertices  and  minimum  vertex 


degree  8  must  have  a  maximum  clique  of  size  at  least 


n  —  8 


and  that  this  bound  is  the 


best  possible  in  tenns  of  n  and  8  . 

Moreover,  we  may  set  a  condition  for  the  existence  of  a  maximum  clique  of  size 
r  in  an  undirected  graph  : 

Observation  1:  Consider  an  arbitrary  graph  G  with  n  vertices.  If  G  contains  a 
clique  of  size  r,  then  there  must  be  at  least  r  vertices,  each  of  which  has  degree  at  least 
r  —  1 . 


Conversely,  if  there  do  not  exist  at  least  r  vertices  of  degree  greater  than  or  equal 
to  r- 1,  then  there  can  be  no  clique  of  order  r  in  the  graph  G.  This  condition  is  used  in  both 
of  our  algorithms  to  set  a  quick  upper  bound  for  the  maximum  clique  size  in  a  given 
graph. 
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2.  An  Upper  Bound  on  the  Number  of  Cliques  of  Size  r 

Let  G  be  an  undirected  graph  with  n  vertices.  Then  an  upper  bound  on  the  largest 
number  of  cliques  of  size  r  (for  some  r  such  that  2  <  r  <  n  )  is  the  number  of  subsets  of 

cardinality  r  of  the  vertex  set  V  ( G ) ,  namely,  f  r  )  .  It  is  obvious  that  this  number  can  be 

Yl 

quite  big  for  very  large  graphs,  especially  if  we  search  for  cliques  of  sizer  withr  «  — .  Of 

course,  besides  using  the  subsets  of  cardinality  r  we  also  have  to  count  the  number  of 
vertices  with  degree  greater  than  or  equal  tor-1,  as  was  discussed  above.  And  so,  in 

(  n\ 

order  for  y,-J  cliques  of  order  r  to  exist,  G  must  have  at  least  r  vertices  of  degree 
greater  than  or  equal  to  r  - 1 .  If  we  denote  the  number  of  nodes  in  G  of  degree  at  least  r 

as  nr,  then  a  better  upper  bound  on  the  number  of  cliques  of  size  r  is  (  /  '  )  . 

B.  GRAPH  THEORY  AND  SOCIAL  NETWORK  ANALYSIS 

As  we  mentioned  above,  a  variety  of  the  tools  and  methods  widely  used  in  this 
thesis  will  come  from  Graph  Theory  and  Operations  Research  areas.  Moreover,  since  the 
applications  we  will  use  come  from  the  Social  Network  Analysis  (SNA)  area  we  will  also 
consider  a  framework  to  introduce  some  basic  concepts  of  the  latter,  and  see  how  all 
these  areas  are  related  to  each  other.  Thus,  we  will  be  able  to  analyze  more  precisely  our 
results,  which  arise  from  the  algorithms’  implementation  on  graphs  which  model  social 
networks,  in  Chapter  V  . 

1.  Graph  Theory 

A  graph  is  an  ordered  pair  consisting  of  two  sets;  the  set  of  vertices,  which 
represent  items  of  interest,  and  the  set  of  edges,  which  connect  any  two  distinct  vertices 
that  satisfy  a  particular  relation  of  interest.  Hence,  a  graph  G  =  (F(G),fi(G))  is 

uniquely  defined  by  its  vertex  set  V(G)  and  its  edge  set  E ( G ) ,  or  by  a  diagram  that 
represents  the  vertices  and  edges  pictorially.  It  can  also  be  described  by  adjacency 
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matrices,  incidence  matrices  or  adjacency  lists,  each  of  which  is  useful  for  particular 
applications  [11,  p  48].  In  our  case,  we  chose  to  describe  our  graphs  by  their  adjacency 
matrices,  which  are  used  as  inputs  to  our  algorithms.  Other  basic  characterizations  of  a 
graph,  or  its  components,  that  we  use  to  solve  the  MCP  are  vertex  degree,  graph  density, 
and  subgraphs,  which  are  derived  by  the  clique  concept  and  its  alternatives.  Below  we 
briefly  introduce  these  graph  theoretic  characteristics  we  used  to  solve  the  MCP. 

a.  Adjacency  Matrix 

As  we  introduced  above,  a  graph  G  =  (V  ( G)  ,E(G))  of  order  n  and  size 

m,  where  F(G)  =  {v1,v2,...,vn}  andis(G)  =  [ex,e2,...,en]  can  be  described  by  its 

adjacency  matrix.  The  adjacency  matrix  of  G  is  the  n  x  n  matrix  A  =  j^a(/ j  ,  where  aif=  1  if 

(v^vj)  is  in  E,  and  is  zero  otherwise.  [1 1].  We  do  not  admit  loops,  and  so  it  follows  that  all 
the  diagonal  elements  of  adjacency  matrix  must  be  zeroes:  au  =  0,  Vi :  1  <  i  <  n  . 

Moreover,  since  we  consider  only  undirected  graphs,  it  follows  that  an 
edge  between  a  vertex  vf  and  another  one  v,  is  also  an  edge  between  v .  and  v,. .  In  other 

words,  aij=aji,\/i  ^  j,  that  is,  the  adjacency  matrix  for  an  undirected  graph  is  a 

symmetric  matrix.  The  adjacency  matrix  of  a  complete  graph  contains  Is  in  all  off- 
diagonal  cells:  a.j  =  ajt  =  1,  Vi  ^  j . 

b.  Vertex  Degree 

The  degree  of  a  vertex  v  in  a  graph  G  is  the  number  of  edges  incident  to  v 
and  is  denoted  by  degG  v ,  or  simply  degv  if  the  graph  G  is  clear  from  the  context  [11]. 

n 

Formally,  the  degree  of  each  vertex  vf  of  a  graph  G  is  degGv(  =  ^a.,  V/  :\<  i  <  n  or 

j= i 

degG  v;  =  |jV(v(.)|,  Vi:  1  <  i<  n  ,  where  N  ( )  denotes  the  set  of  neighbors  of,  or  nodes 
adjacent  to,  a  vertex  v.  e  V (G) . 
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The  degree  of  the  vertices  of  a  graph  G  take  values  that  are  not  arbitrary 
but  are  governed  under  some  basic  rules  and  are  related  to  the  order  and  the  size  of  G. 
Hence,  we  may  introduce  the  First  Theorem  of  Graph  Theory  and  an  observation  about 
the  limits  of  degree  of  the  vertices  as  follows: 

The  First  Theorem  of  Graph  Theory  [11]:  If  G  is  a  graph  of  size  m  , 
then  ^  degG  v  =  2m  . 

veV(G) 

Observation  2:  If  G  is  a  graph  of  order  n  and  v  is  any  vertex  of  G , 
then:0<A(G)<degv<  A(G)<«  -1,  where  8{G )  and  A(G)  denote  the  minimum  and 
the  maximum  degree  of  G,  respectively  [11]. 

c.  Cliques 

We  now  present  the  definitions  related  with  the  clique  and  its  associated 
concepts  accompanied  by  several  examples.  These  concepts  are  strongly  related  with  the 
MCP  and  are  considered  in  different  kinds  of  relaxation  of  the  clique  problem  as  we 
explain  in  the  next  section. 

Given  an  undirected  graph  G(V,E )  with  n  vertices  and  vertex  set  denoted 
by  F(G) ,  a  clique  is  any  complete  subgraph  of  G.  A  clique  of  order  r  is  denoted  by  Kr . 
A  maximal  clique  in  a  graph  G(V,E )  is  a  clique  that  can  not  be  entirely  contained 

within  another  clique  [10],  while  a  maximum  clique  is  the  largest  complete  subgraph  of 
G.  The  order  of  the  maximum  clique  of  G  is  called  the  clique  number  of  G,  and  is 
denoted  by  «(G).  Examples  of  nontrivial  cliques  (a  single  node  or  two  adjacent  nodes 

form  trivial  cliques),  nontrivial  maximal  cliques,  and  the  maximum  clique  of  G  are 
shown  in  Figure  1. 
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Cliques  of  K3 :  {1,  2,  3},  {1,  3,  4}, 

{3.4,6},  {3,4,7}, 
{3,6,7},  {4,  6,7} 
Clique  of  K4  {3,  4,  6,  7} 

Maximal  cliques:  {1, 2,  3}, 
{1,3,4}, 

{3,4,  6,7} 

Maximum  clique  with  co(G)  =  4 : 

{3,4,  6,7} 

Figure  1 .  A  graph  G  and  its  nontrivial  cliques,  maximal  cliques  and  maximum  clique. 

(From  [10]) 

2.  Social  Network  Analysis  (SNA) 

The  pattern  of  relationships  among  the  members  of  a  group  or  a  larger  social 
system  gives  the  relational  structure  among  them.  SNA  is  used  to  study  this  relational 
structure  or  any  other  structural  variables  measured  on  actors  in  the  set  [10].  According  to 
Memon  and  Larsen  in  [19],  “SNA  in  general  studies  the  behavior  of  the  individual  at  the 
micro  level,  the  pattern  of  relationships  (network  structure)  at  the  macro  level,  and  the 
interactions  between  the  two.”  Moreover,  SNA  applies  techniques  to  these  relationships 
and  investigates  how  they  could  be  used  to  infer  more  information  about  the  actors  and 
groups  [20].  More  particularly,  D.  M.  Akbar  presents  the  following  characteristics  that 
SNA  is  intended  to  help  identify  [21]: 

•  Important  individual,  event,  place  or  group. 

•  Dependency  of  individual  nodes. 

•  Leader-Follower  identification. 

•  Bonding  between  nodes. 

•  Vulnerabilities  identification. 

•  Key  players  in  the  network. 

•  Potential  threat  from  the  network. 

•  Efficiency  of  overall  network. 
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In  SNA,  a  subgraph  in  a  graph  is  a  clique  if  it  is  a  maximal  complete  subgraph  of 
three  or  more  vertices.  Mutual  dyads  are  not  considered  to  be  cliques.  For  that  reason  the 
restriction  that  the  clique  contains  at  least  three  nodes  is  included  in  the  definition  [10]. 
We  have  to  notice  that  referring  to  clique  here  Wasserman  and  Faust  mean  maximal 
clique.  Moreover  according  to  Balasundaram  et  al.  in  [5],  “Clique  models  idealize  three 
important  structural  properties  that  are  expected  of  a  cohesive  subgroup,  namely, 
familiarity  (each  vertex  has  many  neighbors  and  only  a  few  strangers  in  the  group), 
reachability  (a  low  diameter,  facilitating  fast  communication  between  the  group 
members)  and  robustness  (high  connectivity,  making  it  difficult  to  destroy  the  group  by 
removing  members).”  However,  the  clique  approach  has  been  criticized  for  its  overly 
restrictive  nature  [5,  10,  22]  and  modeling  disadvantages  [23,  24],  Thus,  alternative 
approaches  have  been  suggested  in  order  to  relax  the  clique  definition  and  different 
models  have  been  developed  in  order  to  relax  different  aspects  of  a  cohesive  subgroup 
[5].  Hence,  we  may  define  the  following  clique  relaxation  models  according  to  the 
property  of  the  cohesive  subgroup  we  want  to  relax  [5,10]: 

1.  Complete  Mutuality  (adjacency):  cliques 

2.  Reachability  and  Diameter  (geodesic  distance):  k-clique 

k-clans 

k-clubs 

3.  Nodal  Degree  (number  of  ties  among  subgroup  members):  k-plexes 

k-cores 

4.  Relatively  Nodal  Degree  (comparison  of  frequency  of  ties  within  to  these 

ones  outside  subgroup):  LS  Sets 

Lambda  Sets 

Thus,  it  is  obvious  that  the  MCP  we  currently  consider  gives  partial 
information  about  the  structure  of  social  networks.  Nevertheless,  it  may  be  further 
considered  under  the  relaxations  presented  above,  in  order  to  provide  a  more  complete 
idea  about  the  network’s  social  structure. 
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III.  SOLVING  THE  MAXIMUM  CLIQUE  PROBLEM 


Once  you  eliminate  the  impossible,  whatever  remains,  no  matter  how 
improbable,  must  be  the  truth. . . . 

Sherlock  Holmes,  by  Sir  Arthur  Conan  Doyle  (1895-1930) 

A.  INTRODUCTION  TO  THE  PRUNING  ALGORITHM 

Our  priming  algorithm  is  based  on  the  “clique  program”  developed  by  Bell  [6], 
which  is  an  algorithm  to  detect  all  maximal  cliques  in  a  graph.  Given  a  graph  G,  the 
pruning  algorithm  uses  Observation  1  to  establish  an  upper  bound  on  the  maximum 
clique  size,  and  then  searches  exhaustively  for  a  clique  of  that  size.  If  no  such  clique  is 
found,  the  upper  bound  is  reduced  by  one,  and  the  algorithm  searches  for  cliques  of  that 
new  size.  If  the  upper  bound  ever  becomes  two,  the  algorithm  halts  (as  each  edge  in  the 
graph  is  a  clique  of  size  two),  otherwise,  it  stops  after  it  has  enumerated  all  cliques  of  the 
maximum  cardinality  in  G. 

Our  pruning  algorithm  is  a  branch-and-bound  algorithm.  As  we  will  explain, 
some  of  the  issues  that  play  a  significant  role  in  the  Pruning  algorithm  are  (1)  the  degree 
of  the  vertices,  (2)  the  number  of  vertices  of  the  same  degree,  (3)  the  adjacency  matrix,  as 
well  as  (4)  its  structure,  its  form  and  its  figure. 

According  to  Observation  1  presented  above  and  considering  the  degree  of  the 
vertices  and  the  number  of  vertices  of  the  same  degree,  the  algorithm  quickly  sets  an 
upper  bound  on  the  clique  number  as  a  first  step.  Given  an  undirected  graph  G{V,E ) 

with  n  vertices,  described  by  its  n  x  n  adjacency  matrix  as  the  input  data,  the  algorithm 
solves  the  MCP  through  the  following  steps  giving  the  corresponding  answers/outputs: 

1 .  Input  the  n  x  n  adjacency  matrix  of  G  , 

2.  Calculate  the  degree  of  each  vertex, 

3.  Set  an  upper  bound  on  the  clique  number,  a>(G ) 


11 


4.  Try  to  find  at  least  one  clique  of  size  co(G) . 

If  such  a  clique  does  not  exist,  reduce  co{G )  and  repeat  this  step  (4). 

Otherwise,  co(G)  is  the  cardinality  of  a  maximum  clique  in  G  . 

5.  Find  all  subsets  of  vertices  that  form  maximum  cliques  for  G  . 

Our  contributions  beyond  Bell’s  [6]  work  is  a  more  efficient  implementation  that 
is  not  limited  to  graphs  of  very  small  cardinality  (Bell’s  code  was  restricted  to  graphs 
with  up  to  six  vertices),  and  the  development  of  the  last  step  of  the  procedure.  There, 
since  we  have  found  that  there  is  at  least  a  maximum  clique  of  a  particular  order,  say 
<z>(G),  we  find  the  maximum  clique(s)  considering  only  those  vertices  with  degree 
greater  than  or  equal  to  co(G)~  1,  and  using  special  data  structures  to  help  avoid 
redundant  or  useless  calculations  in  the  enumeration,  as  presented  in  the  next  section. 

B.  THE  PRUNING  ALGORITHM 

Let  G  =  (F(G),is(G))  be  a  graph  with  «  =  |f(G)|  vertices  and  m  =  \E(G}\ 
edges  and  v(.  e  F(G)  any  vertex  of  G  .  Then  /V  [  v(  ]  is  the  set  of  neighbors  of  vt,  that  is, 
the  set  of  all  adjacent  vertices  to  vt ,  together  with  v;.  itself,  called  the  closed 
neighborhood  of  vt.  The  subgraph  induced  by  N  [  v;  ]  is  denoted  by  (/V  [v;]) .  The 
cardinality  of  /V  [  v;  ]  is  denoted  by  |  /V  [  vv  ] | ;  that  is,  the  number  of  vertices  that  are 

adjacent  to  v;. ,  together  with  v(.  itself.  Similarly,  for  any  v(.  e  V  (G) ,  then  N ( v, )  is  the  set 

of  neighbors  of  v;,  that  is,  the  set  of  all  adjacent  vertices  to  v;. ,  without  v(.  itself,  called 
the  neighborhood  of  vt .  The  subgraph  induced  by  N  (v, )  is  denoted  by  (N  (v(.  .  The 

cardinality  of  /V  ( v;. )  is  denoted  by  N  ( v(  )| ;  that  is,  the  number  of  vertices  that  are 

adjacent  to  v. ,  without  v;  itself. 
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1. 


Algorithm  PRUNING 


1.  Set  the  nxn  adjacency  matrix  A  =  [a;y ]  of  G  ,  where  a,j=\  if  (/,/)  e  E  , 
and  is  zero  otherwise. 

n 

2.  Vv;.  e  V (G)  find  the  degree  of  v. ,  degv,  =  ^r/(/,  V/'  e  [l ,n] ,  and  calculate 

j= i 

the  maximum  degree  in  G,  A(G)  =  max  [deg  v,.] . 

3.  Set  an  upper  bound  on  the  clique  number: 

a.  Set  ®1(G)  =  A(G)  +  1. 

«  u  i 1  '/degv,>®|(G)-l 

b.  Let  inducedSum  =  )f,  where  ^  =(o  otherwise  , 

;=1 

Vi  e  [l,n] . 

c.  If  inducedSum  >  co}  ( G  ) ,  then  go  to  the  next  step  4;  there  might 
be  a  maximum  clique  with  clique  number  co2  (G)  =  cox  (G) , 
where  co2(G)  is  again  a  potential  clique  number  and  forms  an 
upper  bound  for  the  clique  number. 

d.  Else  if  inducedSum  <  co{  (G) ,  then  go  back  to  step  3b  and  set 
col(G)  =  col(G)-l. 

4.  Check  if  there  is  at  least  a  maximum  clique  and  the  clique  number  (o  ( G ) , 
Vv;  e  V (G) :  deg(v!  )  >  co2  (G),  with  i  e  [l ,n] : 

a.  Define  A[v;.] :  A[v;]  =  {jc15jc2,...,jcr}  ,  with  1  <  r  <  |a[v;]|  . 

b.  Create  a  vector,  say  “adjacentVector”,  where  N  [v;]  is 
temporarily  stored. 
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c.  Let  inducedSuml  =  j]  ct  where  ct  =  { {,  if0^^G)  1 ,  Vi  e  [l,  r] . 

i= 1 

d.  If  inducedSuml  >  co2  (G) ,  then  go  to  the  next  step  5;  there  is  at 
least  one  maximum  clique  of  cardinality  <u(G)  <  <x>2  (G) . 

e.  Else  if  inducedSuml  <  co2[G)  then  go  back  to  step  4  and  set 
co2(G)  =  co2(G)-  1. 

5.  Find  the  subset(s)  of  vertices  that  form  the  maximum  clique(s) 
Vv;  e  V (G) :  deg(v;  )  >  <»(G),  with  i  e  [l ,n\ : 

a.  \/vk  e  (  v(. ) ,  with  1  <  A:  <  |  (v.  )|  define  -^(14): 

N (v*  )  =  {f  1,  y 2 , y, } ,  with  1  <  5  <  ( v,  )| . 

b.  V e  iV  ( v; ) ,  with  1  <  k  <  |Ar(vi.  )|  create  a  vector,  say  “stack,” 
where  N  (vk )  is  temporarily  stored. 

c.  For  y,  e  N (vk) ,  that  is  for  the  first  element  of  each  “stack” 
execute  the  following  steps: 

1)  Create  a  vector,  say  “stackMaximumClique,”  where  the 
vertices  that  fonn  a  maximum  clique  will  be  stored,  that  is 
“ stackMaximumClique ”  =  (zx,z2,...,zt} ,  with  1  <  t  <  co(G) . 

2)  Set  Zj  =vi,z2  =vk,z3  =  y, .  These  three  vertices  fonn  a  K3 
clique  so  far. 

3)  Fet  “element A”  =  y,  and  “ element B"  =  y2 .  Check  if 
“element A”  is  adjacent  to  “ element B ”  .  (Do  this  until 
“  element  A”  =  ysl  and  “  element B”  =  ys). 
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a)  If  yes,  then  check  if  “ element B”  is  adjacent  to  each 
element  of  the  “stackMaximumClique”: 


i.  If  yes,  then  add  “ element B ”  in 

“stackMaximumClique”  and  go  back  to  step 
5.C.3),  where  set  “  element  A”  =  y2  and 

“ element B”  =  y3 .  Continue  this  until 

“  element  A”  =  ys]  and  “  element  B"  =  ys . 

ii.  If  no,  go  back  to  step  5.C.3),  where  set 

“elementB”  =  y3 .  Continue  this  until 

“  element  A”  =  ysl  and  “ elementB ”  =  ys . 

b)  If  no,  then  go  back  to  step  5.C.3),  where  set 
“ elementB ”  =  y3 .  Continue  this  until 

“  element  A”  =  ys]  and  “ elementB ”  =  vs . 

Observation  4:  The  largest  size  among  the  “stackMaximumClique”  vectors 
provides  the  clique  number  of  the  graph. 

C.  THE  ENUMERATION  ALGORITHM 

Given  an  undirected  graph  G(V,E )  with  n  vertices,  our  enumeration  algorithm 

finds  all  the  cliques  in  G  through  the  use  of  a  stack  that  always  contains  a  list  of  nodes 
that  comprise  a  clique  of  size  equal  to  the  number  of  nodes  in  the  stack,  top,  and 
successively  attempts  to  add  each  node  not  currently  in  the  stack,  but  adjacent  to  the  top 
node.  If  such  a  node  is  also  adjacent  to  all  other  nodes  in  the  stack,  it  is  added  to  the 
stack,  and  the  algorithm  continues  in  this  manner  until  no  vertices  remain  to  be  added.  It 
then  backtracks,  and  attempts  to  add  the  next  unexplored  node  to  the  remaining  clique. 
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1.  Algorithm  ENUMERATION 


top  =  0 

for  s  =  1  to  n  begin 

STACK[++top]  =  s 
onPath[s]  =  top 
next_node[s]  =  s+1 
while  top  >  0  begin 

i  =  STACK[top] 
while  next_node[i]<=n  begin 
j=next_node[i] 

while  j  <=  n  and  adjacent[i,j]==0 

j++ 

next_node[i]  =j+l 
if  (j  <=  n)  begin 
k=l 

while  k<top  and  adjacent[i,STACK[k]] 
k++ 

if  k==top  begin 

STACK[++top]  =  j 
next_node[j]  =j+l 

end 
else  begin 
top- 

end 

i  =  STACK[top] 

end 

end 

next  s 
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2. 


Discussion 


The  enumeration  algorithm  finds  all  the  cliques  which  exist  in  an  undirected 
graph  and  generates  each  clique  exactly  once.  In  order  to  temporarily  store  the  vertices 
we  have  found  so  far  during  each  turn  and  fonn  a  clique,  we  use  a  one-dimensional  array 
with  (up  to)  n  entries,  STACK\\,  that  is  managed  as  a  stack  using  the  variable  top  as 
pointer  to  the  top  of  the  stack  [25].  Hence,  the  first  top  entries  of  STACK  always  contains 
a  list  of  vertices  that  form  a  clique  of  cardinality  top.  Moreover,  considering  the  vertices 
in  ascending  vertex-number  order,  we  start  investigate  if  each  one  belongs  to  a  clique  by 
setting  it  as  the  first  element  of  the  “stack.”  This  investigation  takes  place  among  this 
particular  vertex  and  those  with  higher  vertex-numbers.  So  the  vertices  which  form  a 
clique  appear  in  ascending  vertex-number  order  on  the  “stack”  and,  thus,  we  avoid 
repetition  of  cliques. 
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IV.  APPLICATION  OF  SOCIAL  NETWORK  ANALYSIS  TO 
CRIMINAL  ACTIVITY:  TERRORIST  NETWORKS 


What  social  network  analysis  contributes  to  counter-terrorism  is  the  ability 
to  map  the  invisible  dynamics  inside  a  terrorist  community. 

P.  V.  Fellman  and  R.  Wright  [17] 


A.  THE  MODEL 

Here  we  develop  a  model  of  a  social  network  to  investigate  our  algorithms’ 
behavior  and  to  examine  and  qualify  their  results.  More  precisely,  we  want  to  model  and 
study  a  real-world  situation.  So,  we  focus  on  terrorist  networks,  which  consist  of  a  direct 
application  area  of  social  networks  [5],  According  to  Balasundaram,  et  ah,  in  [5],  this  “is 
essentially  a  special  application  of  criminal  network  analysis  that  is  intended  to  study 
organized  crimes  such  as  terrorism. ...” 

Throughout  the  analysis  of  a  particular  social  network  and,  more  specifically,  of  a 
terrorist  network,  we  try  to  give  answers  to  crucial  questions  about  its  structure,  like, 
“who  is  (are)  the  leader(s)  and  how  I  could  identify  him(them)?”  or,  “Are  there  any  active 
subgroup(s)  in  the  network  and  how  I  could  recognize  it(them)?”  The  infonnation  that  we 
may  derive  from  a  network  varies,  depending  on  the  type  of  network.  According  to  P.  V. 
Felrnan  et  al.  in  [18],  terrorist  networks  are  “first  and  foremost,  covert,  which  means  that 
they  have  hidden  properties,  and  our  information  about  them  is  incomplete.”  This  will 
become  obvious  in  the  next  section,  where  we  examine  a  real-world  application  of 
terrorist  networks. 

B.  REAL-WORLD  APPLICATION  #  1,  THE  TERRORIST  ATTACK  OF 

SEPTEMBER  11,  2001 

Hence,  as  a  first  application  model  we  selected  a  terrorist  network  of  an  extremely 
tragic  event  which  marked  world  history  and  signaled  the  start  of  the  GWOT  (Global 
War  On  Terrorism).  This  is  the  terrorist  network  which  depicts  the  structure  and  the  links 
among  the  members  who  were  involved  in  the  terrorist  attack  of  September  11,  2001. 
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Valdis  Krebs  used  public  data  that  was  available  before,  but  collected  after  the  event  [5] 
and  constructed  a  graph.  He  presented  this  graph  in  a  remarkable  paper,  [16],  where  he 
made  a  deep  analysis  of  the  terrorist  network  that  caused  this  terrible  attack,  based  on 
SNA.  He  initially  mapped  a  portion  of  the  network  centered  around  the  19  dead  hijackers, 
providing  “some  insight  into  the  terrorist  organization,  yet  it  is  incomplete”  [16].  Krebs 
collected  information  about  the  19  hijackers  and  the  relations  which  connected  them, 
presenting  them  in  a  matrix  named  Early  Hijacker  Matrix,  as  is  shown  below  in  Figure  2 
[16]. 


THE  HIJACKERS  ... 


AND  HOW  THEY  WERE  CONNECTED 


Crashed  into  WTC  (north) 


I1  I  I  □  I  BM 

Crashed  into  WTC  (south) 


A  W 


Mohamed  Atta 

(Egyptian) 

Received  pilot  training 


O 

*-ir 


Marwan  al-Shehhi 
(United  Arab  Emirates) 
Received  pilot  training 


Waleed  M.  Alshehri 

(Saudi) 

Commercial  pilot 


Wail  Alshahri 
(Saudi) 

4  Possible  pilot  training 


n  )  Fayez  Ahmed 

p.ctuff  (Believed  to  be  Saudi) 

«v«itao<e 


Ahmed  Alghamdi 
(Possibly  Saudi) 


Satam  al-Suqami 
(Nationality  unknown) 


no  Abdulaziz  Alomari* 

picture  (Saudi) 

a>*iLatxe  Possible  pilot  training 


Hamza  Alghamdi 

(Believed  to  be  Saudi) 
Possible  pilot  training 


Mohald  Alshehri 
(Nationality  unknown) 
Possible  pilot  training 


Crashed  into  Pentagon 


Khalid  al-Midhar 

(Nationality  unknown) 
Received  pilot  training 

Majed  Moqed 
(Nationality  unknown) 


Salem  Alhamzi* 
(Saudi) 

Possible  pilot  training 


Nawaf  Alhamzi" 
(Saudi) 


Hani  Hanjour 

(Saudi) 


t'j.'.ui.r.iMicg 

Crashed  in  Pennsylvania 

Ziad  Jarrah 
(Lebanese) 

jKyy  Received  pilot  training 

Ahmed  Alhaznawi 
(Saudi) 


Ahmed  Alnami 
(Nationality  unknown) 


Saeed  Alghamdi* 
(Seems  to  be  Saudi) 


’Disputed 

identity 


Attended  same 
technical  college 

Hamburg,  Germany 
Mohamed  Atta 
Marwan  al-Shehhi 
Ziad  Jarrah 


Took  flight  classes 
together 

Pilot  schools 
in  Florida 

Mohamed  Atta 
Marwan  al-Shehhi 

Pilot  schools 
In  San  Diago 

Khalid  al-Midhar 
Nawaf  Alhamzi 

Bought  flight 
tickets  using 
same  address 

•  Mohamed  Atta* 
Marwan  al-Shehhi 
Abdulaziz  Alomari* 

*  Also  used  some 
credit  card 

•  Waleed  M.  Alshehri 
Wail  Alshahri 

•  Fayez  Ahmed 
Mohald  Alshehri 

•  Ahmed  Alghamdi 
Hamza  Alghamdi 


Known  to  be 
together  in  week 
before  attacks 

Stayed  together 
in  a  Florida 
motel 

Mohamed  Atta 
Marwan  al-Shehhi 

Attended  a  gym 
in  Maryland 
(Sept  2-6), 
also  seen  dining 
together 
Khalid  al-Midhar 
Majed  Moqed 
Salem  Alhamzi 
Nawaf  Alhamzi 
Hani  Hanjour 

Bought  flight 
tickets  together 

Mohamed  Atta 
Ziad  Jarrah 
Ahmed  Alhaznawi 

Picked  up  tickets 
bought  earlier  in 
Baltimore 

Khalid  al-Midhar 
Majed  Moqed 

Bought  from  the 
same  travel  agent 
in  Florida 

Ahmed  Alnami 
Saeed  Alghamdi 


Last  known 
address 

Hollywood,  Florida 

Marwan  al-Shehhi 
Waleed  M.  Alshehri 
Wail  Alshahri 
Ziad  Jarrah 
Hani  Hanjour 

Other  cities 
in  Florida 

Mohamed  Atta 
Fayez  Ahmed 
Ahmed  Alghamdi 
Mohald  Alshehri 
Khalid  al-Midhar 
Ahmed  Alhaznawi 
Ahmed  Alnami 
Saeed  Alghamdi 

Outside  Florida 

Satam  al-Suqami 
Hamza  Alghamdi 
Abdulaziz  Alomari 
Majed  Moqed 
Salem  Alhamzi 
Nawaf  Alhamzi 


Figure  2.  Early  Hijacker  Matrix.  (From  [16]) 

After  completing  this  matrix  he  started  building  a  network  considering  more 

actors  and  ties  than  these  ones  among  the  19  hijackers  than  he  had  initially  used.  Finally, 

he  constructed  a  graph  which  mapped  this  terrorist  network,  named  “Hijacker’s  Network 

Neighborhood”  [16]  and  it  is  shown  below  in  Figure  3.  Moreover,  Balasundaram,  et  al., 

presented  this  graph  in  [5],  as  shown  below  in  Figure  4,  having  replaced  the  names  of 
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vertices  by  numbers.  So,  for  practical  reasons,  we  used  this  last  version  of  the  “Hijacker’s 
Network  Neighborhood”  as  the  input  into  our  programs,  since  each  vertex  is  marked  by  a 
number  and  there  is  also  a  correspondence  between  the  names  of  the  hijackers  and  the 
numbers  of  the  vertices.  This  network  corresponds  to  a  graph  with  37  vertices  and  85 
edges,  that  is  a  graph  say  (7(37,85) . 


Wail  Alshehri 


latam  Suqami 


Nabll  al-Marabh 


Waleed  Alshehri 


•ayez  Ahmed 


Mustafa  Ahmed 
Aziz  Al-Omari* 


ihmed  Al  Hazni 


Mam  dan  Darkazanli 


Salem  Aihazmi 


Jarrah 


MounlnEI  Mot 


Khalid 


lamzi  Bin  a] 


lanjour 


iamd{ 

,hmed  Alnami 


med  Alghamdi 


Zacarias  Moussaoi 

w 

Agus  Budiman 

■ 

Ahmed  Khalil  Ibrahim  Sarvjir  Al-Aini 

■  Flight  AA  #11  -  Crashed  into  WTC  North 

■  Flight  AA  #77  -  Crashed  into  Pentagon  Rayed  Motjammad  Abdullah 

■  Flight  UA  #93  -  Crashed  into  Pennsylvania 

■  Flight  AA  #175  -  Crashed  into  WTC  South 

■  Other  Associates  of  Hijackers 


Mohamed  Abdi 


ussattar  Shaikh 


Osama  Awadallah 


Majed  Moqed 


Faisal  Al  Salmi 
Bandar  Aihazmi 


Copyright  €>,  Valdis  Krebs 


Figure  3.  The  “Hijacker’s  Network  Neighborhood.”  (From  [16]) 
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1 


4 


3 


6 


1 

2 

3 

4 

5 

6 

7 

8 
'i 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


Wail  Alshehri 
Satam  Suqami 
Nabil  al-Marabh 
Raed  Hijazi 
Waleed  Alshehri 
Ahmed  Alghamdi 
Mohand  Alsheliri 
Saeed  Alghamdi 
Payez  Ahmed 
Mustafa  Ahmed  Al-Hisawi 
Abdul  Aziz  Al-Omari 
Hamza  Alghamdi 
Ahmed  Alnami 
Ahmed  A1  Haznawi 
Mamoun  Darkazanli 
Mohamed  Abdi 
Marwan  Al-Shehhi 
Zakariya  Essabar 
Salem  Alhazmi 


20  Nawaf  Alhazmi 
21|  Said  Bahaji 

22  Ziad  Jarrah 

23  Mohamed  Atta 

24  Abdussattar  Shaikh 

25  Mounir  El  Motassadeq 

26  Khalid  Al-Mihdhar 

27  Zaearias  Moussaoui 

28  Ramzi  Bin  al-Shibh 

29  Lofti  Raissi 

30  Hani  Hanjour 

31  Osama  Awadallah 

32  Agus  Budiman 

33  Ahmed  Khalil  Ibrahim 
Samir  Al-Ani 

34  Majed  Moqed 

35  Rayed  Mohammed  Abdullah 

36  Faisal  A1  Salmi 

37  Bandar  Alhazmi 


Figure  4.  The  “Hijacker’s  Network  Neighborhood”.  (From  [5]) 


After  running  the  two  programs,  both  of  them  gave  us  the  same  correct  solution  to 
the  MCP.  More  precisely,  they  indicated  that  in  the  “Hijacker’s  Network  Neighborhood,” 
there  is  only  a  maximum  clique  with  clique  number  six  (<»((?)  =  6).  This  maximum 

clique  consists  of  the  following  vertices  17,18,21,22,23,28  and  is  marked  by  a  red  circle 
in  the  graphs  of  Figure  3  and  Figure  4  above.  Moreover,  the  enumeration  algorithm 
provided  all  the  cliques  in  the  graph,  as  shown  on  the  last  pages  of  Appendix  D. 

Analyzing  the  programs’  output  results,  we  may  present  the  following  comments: 

1.  First  of  all  it  is  obvious  that  our  algorithms  run  correctly  and  have 
successfully  implemented  in  the  programming  language  we  selected, 
namely  Java. 

2.  Only  three  out  of  the  six  actors  who  form  the  maximum  clique  were 
among  the  19  hijackers.  The  other  three  actors  where  among  their 
“accomplices  who  did  not  get  on  the  planes”  [16].  The  degrees  of  the  three 
hijackers  of  maximum  cliques  are  the  largest  in  the  whole  network  while 
the  three  accomplices  of  maximum  clique  have  not  so  large  degrees.  The 
six  actors  who  form  the  maximum  clique  in  the  graph  of  “Hijacker’s 
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Network  Neighborhood,”  with  their  corresponding  vertex  numbers  and 
their  degrees,  are  shown  below,  in  Table  1. 


Hijackers 

Accomplices 

Actor’s 

Mohammed 

Marwan 

Ziad 

Zakariya 

Said 

Ramzi 

name 

Atta 

al-Shehhi 

Jarrah 

Essabar 

Bahaji 

Bin 

alShibh 

Vertex 

23 

17 

22 

18 

21 

28 

number 

Degree 

15 

14 

10 

5 

7 

8 

Table  1.  The  actors  who  consist  the  maximum  clique  and  their  degrees 


3.  The  actors  with  the  highest  degree  in  the  graph  (greater  than  or  equal  to 
ten)  are  five,  all  hijackers,  while  three  of  them  are  the  ones  who  belong  to 
the  maximum  clique  and  the  other  two  do  not  belong  to  the  maximum 
clique,  as  is  shown  below,  in  Table  2. 


Hijackers  belonging  to  maximum 

clique 

Hijackers  not  belonging  to 

maximum  clique 

Actor’s 

name 

Mohammed 

Atta 

Marwan 

al-Shehhi 

Ziad  Jarrah 

Nawaf 

Alhamzi 

Hami  Hanjour 

Vertex 

number 

23 

17 

22 

20 

30 

Degree 

15 

14 

10 

10 

10 

Table  2.  The  actors  with  the  highest  degree  (greater  than  or  equal  to  ten) 
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4. 


The  three  hijackers  of  the  maximum  clique  were  pilots  of  a  separate  flight, 
each  one  as  is  shown  below  in  Table  3,  supposing  that  they  had  a 
leadership  role.  In  this  table  they  are  highlighted  in  yellow.  Moreover,  two 
of  the  hijackers  of  the  fourth  flight  (Flight  77,  which  crashed  into 
Pentagon)  do  not  belong  to  the  maximum  clique,  but  they  are  among  the 
actors  with  the  highest  degree  in  the  graph  presented  above.  These  last 
ones  are  highlighted  in  red  in  Table  3.  All  the  others  hijackers  do  not 
belong  to  the  maximum  clique  and  have  intermediate  or  very  small  degree 
in  the  graph. 


American  Airlines 

11  Crashed  into 
WTC  (north) 

United  Aii 
Crashed  in 
(soul 

-lines  11 
to  WTC 
h) 

American  Airlines 
77  Crashed  into 
Pentagon 

United  Airlines  93 
Crashed  into 
Pennsylvania 

Hijacker 

(Vertex 

number) 

Degree 

Hijacker 

(Vertex 

number) 

Degree 

I  lijackel 

(Vertex 

number) 

Degree 

Hijacker 

(Vertex 

number) 

Degree 

Mohamed 

Atta 

(23) 

15 

Marwan 

al-Shehhi 

(17) 

14 

(20) 

10 

Ziad 

Jamah 

(22) 

10 

Waleed  M. 
Alshehri 
(5) 

4 

Hamza 

Alghamdi 

(12) 

6 

(30) 

10 

Saeed 

Alghamdi 

(8) 

6 

Satam  al- 
Suqami 
(2) 

4 

Ahmed 

Alghamdi 

(13) 

3 

Khalid 

al- 

Midhar 

(26) 

4 

Ahmed 

Alhaznawi 

(14) 

3 

Abdulaziz 

Al-Omari 

(ID 

3 

Fayez 

Ahmed 

(9) 

3 

Salem 

Alhamzi 

(19) 

3 

Ahmed 

Alnami 

(13) 

3 

Wail. 

Alshehri 

(1) 

2 

Mohald 

Alshehri 

(7) 

2 

Majed 

Moqed 

(34) 

1 

Table  3.  The  19  Hijackers  and  their  Degrees  in  the  “Hijacker’s  Network 

Neighborhood” 


5.  Mohamed  Atta  seems  to  play  a  significant  role  in  this  terrorist  network 
since  he  belongs  to  the  maximum  clique  and  has  the  highest  degree  of  all 
37  actors  of  the  network. 
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6. 


The  maximum  clique  cannot  give  a  sufficient  and  complete  answer  about 
the  structure  and  the  detection/localization  of  the  leaders  and  the  key 
actors  in  the  network  without  considering  other  parameters  or  metrics  of 
the  graph.  If  we  consider  the  vertices  degree  together  with  the  maximum 
clique,  as  we  analyze  above,  we  may  have  a  deeper  understanding  of  the 
structure  and  leadership  issues  of  the  graph. 

As  strong  evidence  and  confirmation  of  the  above  results,  we  quote  below  some 
comments/conclusions  that  Krebs  drew  in  [16],  where  he  did  a  thorough  analysis  of  the 
“Hijacker’s  Network  Neighborhood”  based  on  SNA  parameters  and  metrics  as  the  main 
centrality  measures,  which  are  Degrees,  Betweenness  and  Closeness: 

•  I  was  amazed  at  how  sparse  the  network  was  and  how  distant  many 
hijackers  on  the  same  team  where  from  each  other. . . 

•  ...in  the  transcript  (Department  of  Defense,  2001)  bin  Laden 
mentions: 

Those  who  were  trained  to  fly  didn’t  know  the 
others.  One  group  of  people  did  not  know  the  other 
group.... 

•  Mohamed  Atta  was  the  ring  leader  of  this  conspiracy. . . 

•  Atta  scores  the  highest  on  Degrees,  and  Closeness  but  not 
Betweenness  centrality.  These  metrics  do  not  necessarily  confirm 
his  leader  status... 

The  above  analysis  is  not  far  from  reality  and,  together  with  the  graph,  it  would  be 
extremely  valuable  if  the  results  were  developed  early  enough  to  prevent  this  tragic  event. 
Unfortunately,  all  the  related  data  was  collected  and  connected  in  a  graph  after  the  attack, 
when  the  analysis  is  easier  since  we  know  what  has  already  happened. 

C.  REAL-WORLD  APPLICATION  #  2,  THE  TERRORIST  ATTACK  OE  U.S. 
EMBASSIES  IN  NAIROBI,  KENYA  AND  DAR  ES  SALAAM,  TANZANIA 
(AUGUST  7, 1998) 

Another  application  model  is  also  a  terrorist  network  of  another  tragic  event 
which  occurred  in  two  countries  of  East  Africa  and  the  targets  were  the  U.S.  Embassies, 
respectively.  On  August  7,  1998,  two  cooperating  Al-Qaeda  cells  carried  out  bombing 
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attacks  against  the  U.S.  Embassies  in  Nairombi,  Kenya  and  Dar-es-Salaam,  Tanzania 
killing  224  people  and  wounding  over  5000  [26],  The  network  which  represents  the  ties 
among  the  members  who  were  involved  in  some  way  in  this  event  is  represented  by  the 
graph  in  Figure  5  [27]  and  is  based  on  data  which  was  collected  after  the  event. 


Figure  5.  The  East  Africa  U.S.  Embassies  Attack  Network  and  the  corresponding 

Maximum  Cliques.  (From  [27]) 

In  this  case  there  is  a  main  difference  to  the  group  structure  comparing  with  this 
one  of  the  9/11  terrorist  attack,  we  presented  above.  Here,  there  are  two  types  of  cells. 
The  first  cell  was  responsible  for  the  preparation  phase  while  the  second  ones  executed 
the  attacks  [26].  This  structure  is  obvious  in  Figure  5,  where  the  subgroups  within  the 
rectangles  denotes  the  two  attack  cells  and  the  rest  actors  form  the  preparation  cell.  This 
network  corresponds  to  a  graph  with  18  vertices  and  51  edges,  that  is  a  graph  say 
G(18,51). 
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After  running  the  two  programs,  both  of  them  gave  us  the  same  correct  solution  to  the 
MCP.  More  precisely  they  gave  that  in  the  “East  Africa  U.S.  Embassies  attack  network” 
there  are  three  maximum  cliques  with  clique  number  five  ( co(G)  =  5 ).  These  maximum 

cliques  are  shown  in  Table  4  and  are  marked  by  a  different  colored  circle  in  the  graphs  of 
Figure  5  above. 


Max  Clique  #1 

Max  Clique  #2 

Max  Clique  #3 

Actor’s  Name 

Degree 

Actor’s  Name 

Degree 

Actor’s  Name 

Degree 

Odeh 

8 

Odeh 

8 

Fahad 

6 

Fazul 

7 

Fazul 

7 

Fadhil 

7 

Al-Owhali 

7 

Al-Owhali 

7 

Awad 

5 

Azzam 

5 

Abdullah 

7 

KKM 

4 

Atwah 

9 

Atwah 

9 

AKG 

4 

Table  4.  The  actors  who  consist  the  maximum  cliques  and  their  degrees 


Analyzing  the  programs’  output  results  we  may  present  the  following  comments  : 

1.  First  of  all  it  is  obvious  that  our  algorithms  run  correctly  and  have 
successfully  implemented  in  the  programming  language  we  selected,  that  is 
Java. 

2.  The  members  of  each  attack  cell  form  or  almost  form  a  maximum  clique.  The 
members  of  the  Nairobi  attack  cell  belong  to  the  first  two  maximum  cliques 
while  the  five  members  of  Dar-es-Salaam  attack  cell  consist  the  third 
maximum  clique.  Hence,  in  this  case,  the  concept  of  maximum  cliques  depict 
extremely  high  valued  cohesive  subgroups  within  the  whole  network, 
assigned  to  execute  the  last  and  most  crucial  act  of  a  terrorist  attack. 

3.  The  number  of  actors  with  the  highest  degree  in  the  graph  (greater  than  or 
equal  to  seven)  is  seven.  In  Table  5  we  may  see  their  role  and  its  position 
within  the  terrorist  network,  as  well  as,  if  they  belong  to  a  maximum  clique  or 
not.  The  information  related  to  the  actor’s  role  and  its  position  within 
network,  have  been  retrieved  by  [26]  and  [27], 
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Actor’s  Name 

Degree 

Actor’s  Role 

Actor’s  Position 

Belonging  to  a 

Maximum 

Clique  or  NOT 

Wahid  el-Hage 

9 

Leader  of  the  East 

African  al-Qaeda  cell 

Preparation  Cell 

NOT 

Matwalli  Atwah 

9 

Fie  acted  as  a 

“communication  bridge” 

among  the  persons  of  the 

two  attack  cells  and  the 

preparation  cell. 

Preparation  Cell 

Max.  Clique 

#1  and  #2 

Mohamed 

Sadeek  Odeh 

8 

Technical  advisor 

responsible  for  carrying 

out  the  bombings. 

Nairobi  Attack 

Cell 

Max.  Clique 

#1  and  #2 

Fazul  Abdullah 

Mohammed 

7 

Planning/  Orchestration/ 

Purchaser 

Fie  was  a  significant 

contributor. 

Nairobi  Attack 

Cell 

Max.  Clique 

#1  and  #2 

Daoud  al- 

Owhali 

7 

Suicide  Bomber. 

Nairobi  Attack 

Cell 

Max.  Clique 

#1  and  #2 

Ahmed 

Abdullah 

7 

Fie  replaced  el-Flage  as 

the  leader  of  al-Qaeda  in 

East  Africa.  Fie  was  the 

“mastermind”  of  the 

coordinated  attack. 

Preparation  Cell 

Max.  Clique 

#2 

Mohammed 

Fadhil 

7 

Operation  Leader. 

(Dar-es-Salaam 

Attack  Cell) 

Max.  Clique 

#3 

Table  5.  The  actors  with  the  highest  degree  (greater  than  or  equal  to  seven) 
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In  the  above  Table  5,  it  is  obvious  that  the  seven  persons  with  the  highest  degree  within 
the  network  played  a  significant  role  to  the  bombing  attacks  and  all  except  one  belong  to 
a  maximum  clique.  The  one  who  does  not  belong  to  a  maximum  clique  is  the  most 
leading  person,  named  Wahid  el  Hage,  having  the  highest  degree  of  all  the  actors.  El- 
Hage  had  “...the  highest  number  of  social  interactions...”  [27]  and  acted  as  the  leader  of 
al-Qaeda  in  East  Africa  [26]  . 

4.  There  were  two  actors  who  belong  to  the  first  two  maximum  cliques, 
without  being  involved  in  the  Nairobi  attack  cell,  which  coincide  with  the 
maximum  cliques  #1  and  #2.  These  ones,  named  Matwalli  Atwah  and 
Ahmed  Abdullah,  were  two  key  persons  who  had  significant  roles  in  this 
attack,  as  it  is  shown  in  Table  5. 

5.  In  this  case,  the  maximum  cliques  give  a  very  good  sense  about  the 
structure  of  this  network.  Together  with  the  number  of  ties  which  is  the 
corresponding  vertex’s  degree,  they  indicate  the  key  persons  who  played  a 
significant  role  in  this  terrorist  attack,  without  considering  any  other  SNA 
metric. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Before  we  fight  the  enemy  we  need  to  'see'  and  better  understand  the 
enemy. ..map  their  networks,  figure  out  their  ’patterns'  for  organizing. 

Valdis  Krebs 


A.  CONCLUSIONS-RESEARCH  SUMMARY 

The  main  question  addressed  in  this  thesis  is: 

Can  we  identify  and  implement  an  exact  algorithm  to  solve  the  Maximum  Clique 

Problem  (MCP)  on  undirected  graphs,  in  a  reasonable  time  frame? 

We  selected  two  algorithms  from  the  current  literature  and  modified  these 
algorithms  to  improve  their  performance,  developing  the  pruning  and  enumeration 
algorithms.  We  then  implemented  both  of  these  in  a  modern,  powerful  and  widely  used 
programming  language,  Java.  Our  testing  involved  applying  both  algorithms  to  real- 
world  situations  which  could  be  modeled  by  undirected  graphs:  terrorist  social  networks. 

We  verified  that  both  algorithms  solve  the  MCP  on  undirected  graphs,  and  are 
quick  on  relatively  small  graphs.  Furthennore,  the  pruning  algorithm  immediately 
establishes  an  upper  bound  on  the  clique  number  of  a  graph  and  then  successively 
improves  this  bound;  thus  the  pruning  algorithm  can  be  tenninated  early  with  a  valid 
upper  bound  on  the  clique  number.  Similarly,  the  enumeration  algorithm  quickly  finds 
small  cliques,  and  successively  discovers  larger  and  larger  cliques  in  the  graph  as  it 
progresses,  each  of  which  provides  a  lower  bound  on  the  clique  number  of  the  graph.  It 
can  be  terminated  early  with  a  valid  lower  bound  on  the  clique  number. 

If  both  are  run  simultaneously,  each  provides  a  bound  the  other  cannot,  and  an 
interval  of  uncertainty  can  be  established  that  will  eventually  be  reduced  to  zero,  at  which 
point  a  maximum  clique  will  have  been  found. 

Moreover  analyzing  a  social  network  using  the  SNA  methods  and  measures,  as 
we  extensively  presented  in  Chapters  I  and  IV,  maximum  cliques  and  vertices  degree 
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may  provide  enough  information  about  the  social  structure  of  the  network  under 
investigation.  More  particularly,  the  maximum  clique(s)  and  the  clique(s)  as  well,  may 
give  us  the  most  cohesive  subgroups,  while  the  vertices  with  the  largest  vertex-degrees 
may  show  the  most  central  actors  in  a  social  network.  However,  since  cliques  have  been 
criticized  for  their  restrictive  nature  as  we  precisely  explain  in  Chapter  I,  we  may 
consider  and  study  one  of  the  branches  of  the  relaxation  clique  problem  and  analyze  a 
social  network  under  this  concept.  This  may  be  the  subject  for  an  extension  of  this  current 
thesis  and  an  area  for  further  investigation.  It  is  worthwhile  to  point  out  that  each 
algorithm  can  be  modified  in  fairly  straightforward  ways  to  allow  various  relaxations  of 
the  definition  of  a  clique. 

B.  FUTURE  RESEARCH 

After  presenting  the  conclusions  above,  it  is  obvious  that  there  are  many  areas  and 
aspects  of  this  problem’s  approach  for  improving  with  further  study  and  future  research. 
Some  suggestions  are  provided  below: 

1.  Since  the  pruning  algorithm  generates  the  maximal  cliques  with 
repetitions,  it  would  be  more  efficient  to  be  modified  in  such  a  way  to 
eliminate  them  and  generate  the  cliques  or  the  maximal  cliques  or  the 
maximum  cliques  exactly  once. 

2.  Both  computing  codes  consider  as  the  input  data  the  under  investigation 
graph’s  adjacency  matrix.  This  method  is  not  so  practical  for  very  large 
graphs,  where  the  user  have  to  set  thousands  of  entries.  Hence  a 
modification  to  the  computing  code  in  order  to  read  and  set  the  adjacency 
matrix  from  a  file  text  would  make  it  more  practical  and  efficient. 

3.  Both  algorithms’  concept  to  solve  the  MCP  is  based  on  the  adjacency 
matrix  of  the  under  investigation  graph.  Someone  may  consider  another 
type  of  data  structure  to  describe  the  graph  and  make  the  algorithms  run 
faster.  Such  data  structure  may  be  as  follows: 
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a.  To  use  a  forward-star  data  structure  in  addition  to  the  adjacency 
matrix  in  order  to  make  finding  the  next  adjacent  vertex  faster. 

b.  To  use  a  forward-star  data  structure  where  the  adjacency  list  of 
each  vertex  would  be  guaranteed  to  be  sorted  in  ascending  order  of 
tail  vertex  number. 

Since  cliques  have  been  criticized  for  their  restrictive  nature  as  we 
precisely  explain  in  Chapter  II,  someone  may  consider  and  study  one  of 
the  branches  of  the  relaxation  clique  problem,  extend  the  current 
algorithms  and  codes  to  this  direction  and  analyze  a  social  network  under 
this  concept. 
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APPENDIX  A.  NOTATIONS 


G  =  {V{G),E{G )) 

V, V(G) 

E,E(G ) 
n 

m 

degG  v ,  degv 
M(v) 

(M(v)) 

JvM 

GM) 

H^G 

(s),(x) 

W ,  W(u,v ) 

C„ 

Kn 

dG  (w,v) ,  d  ( u,v ) 
diam  ( G ) 


graph  consisting  of  n  =  F (G)|  vertices  and  m  =  /:  (G’)|  edges. 

vertex  set  of  a  graph  G  . 

edge  set  of  a  graph  G  . 

order  of  G  ,  n  =  \v  (G)| 

size  of  G  ,  m  =  |fs  (G)| 

degree  of  a  vertex  v  in  a  graph  G  . 

neighborhood  of  a  vertex  v ,  | N  ( v)|  =  deg  v 

subgraph  induced  by  N  (v) 

closed  neighborhood  of  a  vertex  v 

subgraph  induced  by  [v] 

II  is  a  subgraph  of  a  graph  G  ,  or  G  contains  H  as  a  subgraph, 
induced  or  edge-induced  subgraph 
u  —v  walk  in  G 
path  of  order  n 
cycle  of  order/; 
complete  graph  of  order  n 
distance  between  u  and  v 
diameter  of  G 
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G 

complement  of  a  graph  G 

Kn 

empty  graph  of  order  n 

S(G) 

minimum  degree  of  G 

A  (G) 

maximum  degree  of  G 

A  =  [a.._ 

j  adjacency  matrix  of  G 

B=W 

incidence  matrix  of  G 

cd(G) 

clique  number  of  G 

w(e) 

weight  of  an  edge  e  of  G 

g,  =  g2 

G,  is  isomorphic  to  G2 

D 

density  of  G 

CD  w 

degree  centrality  of  a  vertex  v  of  G 

QW 

closeness  centrality  of  a  vertex  v  of  G 

c.W 

betweenness  centrality  of  a  vertex  v  of  G 
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APPENDIX  B.  DEFINITIONS  RELATED  TO  GRAPH  THEORY 


The  following  terms  provide  a  core  working  vocabulary  for  discussing  Graph 
Theory.  The  definitions  for  these  terms  are  derived  from  rill. 

Adjacency  matrix:  Let  Gbe  a  graph  of  order  n  and  size  m,  where 
F(G)  =  {v15v,,...,vn}  andfs(G)  =  [el,e2,...,en] .  The  adjacency  matrix  of  G  is  the 

n  x  n  matrix  A  =  [atj  ]  ,  where  a..  =  { '0  • 

Adjacent  vertices  inG  :  two  vertices  are  called  adjacent  in  G  if  there  is  an  edge 
between  them.  In  other  words  if  the  pair  of  these  vertices  is  an  element  of  the  edge 
setE^G ) . 

Bridge:  an  edge  e  =  uv  of  a  connected  graph  G  is  called  a  bridge  of  G  if  G-e 
is  disconnected. 

Circuit:  a  circuit  in  a  graph  is  a  closed  trail  of  length  3  or  more. 

Clique:  a  clique  in  a  graph  G  is  a  complete  subgraph  of  G  . 

Clique  number:  the  order  of  the  largest  clique  in  a  graph  G  and  it  is  denoted 
by  <»(G) . 

Closed  neighborhood:  the  set  of  neighbors  of  a  vertex  v  together  with  the  vertex 
v  itself  and  is  denoted  by  N[v] .  In  other  words,  N[v]  =N  (v)  U  {v}  . 

Closed  walk:  a  u  -  v  walk  W  in  G  where  u  =  v . 

Complement:  the  complement  of  a  graph  G  ,  denoted  by  G  ,  is  that  graph  whose 

vertex  set  is  V  ( G )  and  such  that  for  each  pair  u,  v  of  vertices  ofG  ,  uv  is  an  edge  of  G  if 
and  only  if  uv  is  not  an  edge  of  G  .  Observe  that  if  G  is  a  graph  of  order  n  and  size  m  , 

then  G  is  a  graph  of  order  n  and  size  ( n2 )  -  m  . 
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Complete:  a  graph  G  is  complete  if  every  two  distinct  vertices  of  G  are  adjacent. 
A  complete  graph  of  order  n  is  denoted  by  Kn .  Therefore,  Kn  has  the  maximum  possible 
size  for  a  graph  with  n  vertices. 

Connected:  let  u  and  v  be  two  vertices  of  a  graph  G  .  Then  u  is  connected  to  v 
(and  also  v  is  connected  to  u)  if  G  contains  a  u-v  path  in  G  .  So,  saying  that  u  and  v 
are  connected  only  means  that  there  is  some  u-v  path  in  G  ;  it  doesn’t  say  that  u  and  v 
are  necessarily  joint  by  an  edge. 

Connected  graph:  a  graph  G  is  connected  if  every  two  vertices  of  G  are 
connected. 

Cut-vertex:  a  vertex  vin  a  connected  graph  G  is  a  cut-vertex  of  G  if  G-v  is 
disconnected. 

Cycle:  is  a  circuit  that  repeats  no  vertex,  except  for  the  first  and  last.  A  k-cycle  is 
a  cycle  of  length  k. 

Degree:  the  degree  of  a  vertex  v  in  G  is  the  number  of  edges  incident  to  v  and  is 
denoted  by  degG  v ,  or  simply  deg  v  if  the  graph  G  is  clear  from  the  context. 

Degree  sequence:  if  the  degrees  of  vertices  of  a  graph  G  are  listed  in  a  sequence 
5,  then  5  is  called  a  degree  sequence  of  G  . 

Density  of  a  graph:  it  is  the  proportion  of  possible  edges  that  are  actually  present 
in  the  graph.  It  is  the  ratio  of  number  of  edges  present  to  the  maximum  possible  and  is 
denoted  by  D . 

Diameter:  the  greatest  distance  between  any  two  vertices  of  a  connected  graph 
G  is  called  the  diameter  of  G  and  is  denoted  by  diam  ( G) . 

Digraph  (or  Directed  graph):  is  a  finite  nonempty  set  V  of  objects  called 
vertices  together  with  a  set  E  of  ordered  pairs  of  distinct  vertices. 

Disconnected  graph:  a  graph  G  is  called  disconnected  if  it  is  not  connected. 


38 


Distance:  the  distance  between  u  and  v is  the  smallest  length  of  a  u-v path  in 
G  and  is  denoted  by  dG  (w,v)  or  simply  by  d (w,v)  if  the  graphG  under  consideration  is 

clear.  Hence  if  d(u,v)  =  k,  then  there  exists  a  u  -  v path P  :u  =  v0,vl,...,vk  =v,  of  length 
k  in  G  ,  but  no  u-v  path  of  smaller  length  exists  in  G . 

Edge:  is  a  2-element  subset  of  the  vertex  set  L(G) .  Edges  are  sometimes  called 

lines. 

Edge-induced  subgraph  ofG:  is  the  subgraph  (x)  induced  by  X  which  has 
edge  set  X  and  consists  of  all  vertices  that  are  incident  with  at  least  one  edge  in  X 

Edge  set:  is  the  set  E  (G)  consisting  of  all  the  edges  of  a  graph  G. 

Empty  graph:  the  graph  that  has  n  vertices  and  no  edges  is  called  the  empty 
graph  of  order  n  and  is  denoted  by  Kn  . 

Geodesic  path:  a  u-vpath  of  length  d(u,v)  is  called  a  u-v  geodesic. 

Graph:  a  graph  G  =  (V(G),fs(G))  is  an  ordered  pair  of  two  sets  L(G)  and 
E  ( G ) ,  where  V  is  a  finite  nonempty  set  of  objects  called  vertices  (the  singular  is  vertex ) 
and  E  is  a  set  of  2-element  subsets  of  V  called  edges.  One  could  also  use  G  =  (V,E).  At 
times,  it  is  useful  to  write  L  ( G )  and  E  ( G )  rather  than  V  and  E  to  emphasize  that  these 
are  the  vertex  and  edge  sets  of  a  particular  graphG  .  The  graph  G  =  (V (G),E(G)J  is 
consisting  of  n  =  L  (G)|  vertices  and  m  =  G  (G)|  edges. 

Incidence  matrix:  Let  G  be  a  graph  of  order  n  and  size  m,  where 
F(G)  =  {v15v2,...,vn}  and  E(G)  =  {el,e2,...,em}  .  The  incidence  matrix  of  G  is  the  nxm 

,  •  n  T  j  1  t  r  ( 1  if  v,  is  incident  with  e; 

matrix  B  =  \_by  J  ,  where  b..  =  { 0  otherwise 
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Induced  subgraph:  a  subgraph  F  of  a  graph  G  is  called  an  induced  subgraph  of 
G  if  whenever  u  and  v  are  vertices  of  F  and  uv  is  an  edge  of  G  ,  then  uv  is  an  edge  of 
F  as  well. 

Isomorphism:  it  is  the  one-to-one  correspondence  (f>  from  V (G,)  to  V(G2 )  such 
that  u  jVj  e  E  (Gt )  if  and  only  if  cp(ux )  (p{vx )  e  E  (G2 ) . 

Isomorphic  graphs:  two  graphs  Gx  and  G2  are  isomorphic  (have  the  same 
structure)  if  there  exists  a  one-to-one  correspondence  (j)  from  V  ( G] )  to  V  ( G2 )  such  that 
MjVj  e  E ( G] )  if  and  only  if  (p{ux ) <7? ( j  e  E(G2).  If  Gx  and  G2  are  isomorphic  graphs, 
then  we  say  that  G,  is  isomorphic  to  G2  and  we  write  Gx  =  G2 . 

Length  of  a  path:  is  the  number  of  edges  encountered  in  a  path. 

Length  of  a  walk:  is  the  number  of  edges  encountered  in  a  walk  (including 
multiple  occurrences  of  an  edge  if  used  in  the  walk). 

Maximal:  a  subgraph  is  said  to  be  maximal  with  respect  to  some  property  if  that 
property  holds  for  the  subgraph,  but  does  not  hold  if  additional  vertices  and  the  edges 
incident  with  them  are  added  to  the  subgraph  [10]. 

Maximal  clique:  a  maximal  clique  in  a  graph  G  is  a  clique  that  can  not  be 
entirely  contained  within  another  clique  [10]. 

Maximum  clique:  a  maximum  clique  in  a  graph  G  is  the  largest  complete 
subgraph  of  G . 

Maximum  degree:  the  maximum  degree  ofG  is  the  maximum  degree  among  the 
vertices  ofG  .  It  is  denoted  by  A(G) . 

Minimum  degree:  the  minimum  degree  of  G  is  the  minimum  degree  among  the 
vertices  of  G  .  It  is  denoted  by  S(G) . 

Multigraph:  consists  of  a  finite  nonempty  set  V  of  vertices  and  a  set  E  of  edges, 
where  every  two  vertices  are  joined  by  a  finite  number  of  edges  (possibly  zero) 
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Neighborhood  of  v:  is  the  set  7V(v)  of  neighbors  of  a  vertex  v.  The 
neighborhood  of  v  is  sometimes  called  open  neighborhood  of  v.  The  cardinality  of  the 
neighborhood  of  v  equals  to  the  degree  of  v,  that  is  deg  v  =  TV  (v)| . 

Neighbors:  two  vertices  are  called  neighbors,  if  they  are  adjacent  in  G. 

Non-isomorphic  graphs:  if  two  graphs  G,  and  G2  are  not  isomorphic,  then  they 
are  called  non-isomorphic  graphs. 

Non-regular  graph:  a  graph  G , which  at  least  two  vertices  have  not  the  same 
degree,  is  called  non-regular. 

Nontrivial  graph:  is  a  graph  with  order  at  least  2. 

Open  walk:  is  a  u-v  walk  Win  G  where u  ^  v  . 

Order  ofG :  is  the  number  of  vertices  in  G  ,  in  other  words  the  cardinality  of  the 
vertex  setV (G) .  That  is  the  order  of  G  is  n  =  \v (G)|  . 

Path:  a  u-v  path  Win  G  is  a  w-vwalk  in  which  no  vertices  is  repeated.  If  no 
vertex  in  a  walk  is  repeated,  then  no  edge  is  repeated  either.  Hence  every  path  is  a  trail, 
while  not  every  trail  is  a  path. 

Proper  subgraph:  a  graph  H  is  called  a  proper  subgraph  of  a  graph  G  if 
flcG  and  either  V  (//)  cz  V  (G)  or  cz  fs(G) . 

Regular:  a  graph  G  ,  where  S  ( G)  =  A  ( Gj  ,  that  is  all  the  vertices  of  G  have  the 
same  degree,  is  called  regular.  If  deg  v  =  r  for  every  vertex  v  of  G,  where  0  <r<  n  - 1 , 
then  G  is  r-regular  or  regular  of  degree  r. 

Size  of  G  :  is  the  number  of  edges  in  G  ,  or  the  cardinality  of  the  edge  setV(G). 
That  is  the  size  of  G  is  m  =  |is(G)| . 

Spanning  subgraph:  if  a  subgraph  H  of  a  graph  G  has  the  same  vertex  set  asG  , 
that  is  V (//)  =  V (Gj ,  then  II  is  a  spanning  subgraph  ofG  . 
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Subgraph:  a  graph  H  is  called  a  subgraph  of  a  graph  G  ,  written // c  G  ,  if 
V (//)  cz  V (G)  and  £(//)  c:  2?(G) .  We  also  say  that  G  contains  H  as  a  subgraph. 

Subgraph  of  G  induced  by  S  :  if  5*  is  a  nonempty  set  of  vertices  of  a  graph  G  , 
then  the  subgraph  of  G  induced  by  S  is  the  induced  subgraph  with  vertex  set  S .  It  is 
denoted  by  (S')  or  by  (S')  to  emphasize  that  this  is  an  induced  subgraph  of  G 

Trail:  a  u  —  v  trail  W in  G  is  a  u-v  walk  in  which  no  edge  is  traversed  more  than 
once,  that  is  no  edge  is  repeated,  while  no  such  condition  is  placed  on  vertices. 

Trivial  graph:  is  a  graph  with  exactly  one  vertex. 

Trivial  walk:  is  a  walk  of  length  0. 

Unweighted  graph:  a  graph  G  each  of  whose  edges  is  not  assigned  a  number 
(called  cost  or  weight  of  the  edge)  forms  an  unweighted  graph. 

Vertex  (Vertices):  is  a  combinatorial  element  in  terms  of  which  a  graph  is 
defined.  Vertices  are  sometimes  called  points  or  nodes. 

Vertex  set:  is  the  set  V  (G)  consisting  of  all  the  vertices  of  a  graph  G  . 

Walk:  a  u  -v  walk  W  in  G  is  a  sequence  of  vertices  in  G  ,  beginning  with  u  and 
ending  at  vsuch  that  consecutive  vertices  in  the  sequence  are  adjacent.  That  is  to  say  that 
W  can  be  expressed  as  W:  u  =  v0,vl,...,vk  =  v ,  where  k  >  0  and  v,.  and  vM  are  adjacent  for 

i  =  0, 1, 2, ..., k  - 1 .  Each  vertex  v;.  (0  <  z  <  k)  and  each  edge  vtvM  (0  <  i  <  k  - 1)  is  said  to 
lie  on,  or  belong  to  W . 

Weighted  graph:  a  graph  G  each  of  whose  edges  is  assigned  a  number  (called 
cost  or  weight  of  the  edge)  forms  a  weighted  graph.  The  weight  of  an  edge  e  of  G  is 
denoted  by  w{e) . 
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APPENDIX  C.  SOCIAL  NETWORK  ANALYSIS  DEFINITIONS 


The  following  terms  provide  a  core  working  vocabulary  for  discussing  Social 
Network  Analysis.  The  definitions  for  these  terms  are  derived  by  f  1 01 . 


Actor:  it  is  the  social  entity.  Actors  are  discrete  individual,  corporate,  or 
collective  social  units. 

Actors  set:  it  is  the  entire  collection  of  actors  on  which  we  take  measurements. 

Affiliation  Network:  it  is  a  two-mode  social  network,  in  which  only  one  set  has 
actors,  while  the  second  mode  is  a  set  of  events  to  which  the  actors  belongs. 

Centrality:  the  centrality  determines  the  importance  of  an  actor  within  the 
network.  It  is  denoted  by  CA  (vj ,  where  v(.  eV(G)  and  A  is  a  generic  measure. 

Cohesive  subgroup:  a  cohesive  subgroup  is  a  subset  of  actors  among  whom  there 
are  relatively  strong,  direct,  intense,  frequent  or  positive  ties. 

Dichotomous  relations:  they  are  those  binary  relations  that  are  coded  as  either 
present  or  absent  for  each  pair  of  actors. 

Directional  relation:  it  is  the  relation  where  the  relational  tie  between  a  pair  of 
actors  has  an  origin  and  a  destination;  that  is  the  tie  is  directed  from  one  actor  in  a  pair  to 
the  other  actor  in  the  pair. 

Dyad:  a  dyad  consists  of  a  pair  of  actors  and  the  (possible)  tie(s)  between  them.  It 
is  the  unit  of  social  network  analysis. 

Event:  the  events  are  often  defined  on  the  basis  of  membership  in  clubs  or 
voluntary  organizations,  attendance  at  social  events,  sitting  on  a  board  of  directors  or 
socializing  in  a  small  group. 

Group:  a  group  is  a  collection  of  all  actors  on  which  ties  are  to  be  measured. 
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Mode:  a  mode  is  a  distinct  set  of  entities  on  which  the  structural  variables  are 
measured. 

Multiple  Relations:  they  are  the  case  where  more  than  one  relations  are 
measured  on  a  single  set  of  actors. 

Non-directional  relation:  it  is  the  tie  between  a  pair  of  actors  which  does  not 
have  a  direction. 

One-mode  network:  it  is  a  social  network  where  all  actors  are  all  of  the  same 
type,  i.e.  people  in  a  work  group.  It  is  the  most  common  type  of  network,  since  all  actors 
come  from  one  set. 

Prestige:  it  is  a  property  that  characterizes  an  actor  in  a  directed  network  with 
directional  relations  and  qualify  an  actor  by  the  number  of  indegrees  ties  rather  than  the 
outdegrees  ones. 

Relation:  it  is  the  collection  of  ties  of  specific  kind  among  members  of  a  group. 
Relational  tie:  a  tie  establishes  a  linkage  between  a  pair  of  actors.  It  is  also  called  social 
tie. 

Single  Relation:  it  is  a  relation  where  each  actor  in  the  actor  set  relates  to  every 
other  actor  of  this  relation. 

Social  Network:  a  social  network  consists  of  a  finite  set(s)  of  actors  and  the 
relations  defined  on  them. 

Social  Network  data:  consists  of  one  (or  more)  relations  measured  among  a  set 
of  actors. 

Subgroup:  a  subgroup  of  actors  is  any  subset  of  actors  and  all  (possible)  ties 
among  them. 

Triad:  it  is  a  subset  of  three  actors  and  the  (possible)  tie(s)  among  them. 

Valued  relations:  they  are  those  relations  which  can  take  a  range  of  values, 
indicating  the  strength,  intensity,  or  frequency  of  the  ties  between  each  pair  of  actors. 
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APPENDIX  D.  THE  CODE  OF  THE  ENUMERATION  ALGORITHM 


public  class  All_Cliques_Enum_Adj_Terror_Thesis  { 

/* 

*  File:  Maximum  Clique 

* 

*  Created  on  April  03,  2008 

* 

*  1st  Modified  on  April  06,  2008 

*  2nd  Modified  on  April  28,  2008 

* 

*  *************THIS  WORKS  FOR  the  Terror  Graph  G(37,85)  WITH  A  K6 
MAX  CLIQUE,  BASED  ON  ENUMARAT I ON  METHOD  ***** 

* 

* 

*/ 

public  static  void  main ( String [ ]  args)  { 

//  initialize  instance  variables 
int  vectorsNumber  =  37; 

int [ ]  degree  =  new  int [vectorsNumber ] ; 

int [ ]  booleanDegree  =  new  int [vectorsNumber ] ; 

int  maxCliqueDegree  =  0; 

int  tempMaxCliqueDegree  =  0;//  This  is  equal  to  (Potential 
Maximum  Clique  #  -1),  (PMC#-1) .  Its  value  decreases  upon  the  condition 

which  are  met,  until  Maximum  Clique  is  found.  Then  Clique#= (updated  by 
the  algorithm) PMC# 

int  initialTempMaxCliqueDegree  =  0;//  This  is  equal  to  the 
initial  value  of  the  (Potential  Maximum  Clique  #  -1),  (PMC#-1) .  Its 

value  remains  constant. 

int  newTempMaxCliqueDegree  =  0; 
boolean  maxClique  =  false; 
boolean  initialMaxClique  =  false; 

int  cliqueNumber  =  -1;//  this  is  the  clique  number 
int[][]  table  =  new  int [vectorsNumber ] [vectorsNumber ]; // 

This  is  the  nxn  adjacency  matrix  of  the  given  graph  G  (with  n  vertices) . 

int [ ]  count  =  new  int [vectorsNumber ]; //  counter  how  many 
Maximum  Cliques  there  exist 

int [ ] []  maximumClique  =  new  int [40] [vectorsNumber];//  This 
int[] [] []  clique  =  new 
int [vectorsNumber ] [1140] [vectorsNumber] ; 
int  row,  column; 

// 

************************************************************************ 

**************************************************** 

//  I.  PUT  AN  UPPER  BOUND 

//  1.  Set  the  nxn  adjacency  matrix  of  the  given  graph  G 
(with  n  vertices) . 

table  [0]  [0]  =  0; 
table  [0]  [1]  =  1; 
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table  [0]  [2]  =  0; 
table  [0]  [3]  =  0; 
table  [0]  [4]  =  1; 


table [36] [33]  =  0; 
table [36] [34]  =  1; 
table  [36]  [35]  =  0; 
table [36] [36]  =  0; 


System. out. println ("Maximum  Clique  Problem  of  G("  + 
vectorsNumber  +  ",85)"); 

System. out . println ( )  ; 

//  2.  Find  the  degree  of  each  vertex  by  adding  the  entries 
of  each  row  of  the  adjacency  matrix. 

System. out . println ( "Find  the  degree  of  each  vertex  "); 
for  (row  =  0;  row  <  vectorsNumber;  row++)  { 

int  sum  =  0;//  This  is  the  degree  of  each  vertex 
for  (column  =  0;  column  <  vectorsNumber;  column++) 
sum  =  sum  +  table [row] [column] ; 
degree [row]  =  sum; 

} 

/ /  print  the  degree  of  each  vertex 

for  (row  =  0;  row  <  vectorsNumber;  row++)  { 

System. out . println ( "degree (v"  +  (row  +1)  +  "  )=  "  + 

degree [row] ) ; 

} 

System. out. print ln() ; 

//  3.  Find  the  Potential  Max  Clique#,  say  PMQ#, 

//  (Potential  Max  Clique#  =  max  Vertex  Degree+1) . 
for  (row  =  0;  row  <  vectorsNumber;  row++) 

if  (degree [row]  >  tempMaxCliqueDegree)  { 
tempMaxCliqueDegree  =  degree [row]; 

} 

System. out . println (" Potential  Max  Clique  Degree:  " 

+  (tempMaxCliqueDegree  +  1)  + 

initialTempMaxCliqueDegree  =  tempMaxCliqueDegree; 


//  4.  Let  v  k  the  vertices  with  degree  greater  than  or  equal 
to  ( PMC#-1 ) , f or  some  k,  l<=k<=n.  Find  the  #  of  vertices  v  k  with  degree 
greater  than  or  equal  to  (PMC#-1),  say  "inducedSum" . 

do  { 

int  inducedSum  =  0;//  the  #  of  vertices  with  degree 
greater  than  or  equal  to  (PMC#-1) 

for  (row  =  0;  row  <  vectorsNumber;  row++)  { 

if  (degree [row]  >=  tempMaxCliqueDegree)  { 
booleanDegree [row]  =  1; 
inducedSum  =  inducedSum  + 

booleanDegree [row] ; 

}  else 
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} 


booleanDegree [row] 


0; 


//  4a.  If  "inducedSum"  is  greater  than  or  equal 
PMC#,  then  maybe  there  is  a  Max  Clique  with  Clique#  =  PMC#. 

if  (inducedSum  >=  tempMaxCliqueDegree  +1)  { 

System. out .println ( "Maybe  there  is  a  K" 

+  (tempMaxCliqueDegree  +1)  + 


Clique . " )  ; 


System. out .println  ( )  ; 
initialMaxClique  =  true; 


to 


"  Max 


//  4b.  Else  if  "inducedSum"  is  NOT  greater  than 
or  equal  to  PMC#,  there  is  not  a  Maximum  Clique  with  Clique#  =  PMQ# .  Set 
PMQ#=PMQ#-1  and  continue  by  going  back  to  step  4. 

}  else  { 


false) ; 


System. out . println ( "There  is  not  a  K" 

+  (tempMaxCliqueDegree  +1)  +  "  Max  Clique."); 
System. out . println ( "Check  for  K" 

+  (tempMaxCliqueDegree)  +  "  Max  Clique."); 
System. out. print ln() ; 
tempMaxCliqueDegree-- ; 

} 

}  while  (tempMaxCliqueDegree  >  0  &&  initialMaxClique  == 

// 


************************************************************************ 


********************************************* 

int  top  =  -1; 

for  (row  =  0;  row  <  vectorsNumber ;  row++)  { 


int [ ]  stack  =  new  int [vectorsNumber] ; //  This 
//  is  the  vector  where  the  vertices  of  the 
//  neighborhood  of  each  element  of  N[V_k]  are 
//  stored. 


int [ ]  nextNode  =  new  int [vectorsNumber ];/ /  This 
//  is  the  vector  where  the  vertices  consisting  a 
//  Maximum  Clique  are  stored. 

stack [++top]  =  row; 
nextNode [row]  =  row  +  1; 


while  (top  >=  0)  { 

int  i  =  stack[top]; 

while  (nextNode [i]  <=  vectorsNumber)  { 
int  j  =  nextNode [i] ; 

while  (j  <  vectorsNumber  &&  table [i] [j]  == 

0)  { 

j++; 

} 

nextNode [i]  •-  j  11; 
if  (j  <  vectorsNumber)  { 
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==  1)  { 


tempMaxCliqueDegree;  spyros++) 


int  k  =  0 ; 

while  (k  <  top  &&  table[j] [stack[k]] 
k++ ; 

} 

if  (k  ==  top)  { 
top++; 

stack [top]  =  j; 
nextNode[j]  =  j  +  1; 

System. out . print In ( ) ; 

System. out .print In ( ) ; 

if  (top  >  cliqueNumber )  { 

cliqueNumber  =  top; 

} 

for  (int  spyros  =  1;  spyros  <= 
if  (top  ==  spyros)  { 


System. out .println ( ) ; 

System. out. println ("There  is  a  K"  +  (top  +1)  +  "  Clique."); 


System. out . println ( "The  Clique  is:  "); 

for  (int  m  =  0;  m  <=  top;  m++)  { 


System. out .print ( (stack [m]  +1)  +  "  "); 

clique [top  +  1] [count [top  +  1]] [m]  =  stack[m] ; 

} 


System. out . println  ( ) ; 


count [top  +  1]++; 


+  1 )  )  ; 


} 


} 

} 

}  else  { 

top--; 

} 

if  (j  <  vectorsNumber )  { 

i  =  stack [top] ; 

} 


} 

int [ ]  maximumCliquel  =  new  int [cliqueNumber  +  1]; 

//  print  the  outputs 
System. out. println() ; 

System. out. println ("The  Maximum  Clique  is  K"  +  (cliqueNumber 


System. out. println ("There  is  "  +  count [cliqueNumber  +  1] 
+  "  Maximum  Clique"); 
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+  1)  +  " 
1]  [n]  [m] 


System. out. println ("The  Maximum  Clique  is:  "); 

//  for  (int  k  =  0;  k  kvectorsNumber ;  k++)  { 

for  (int  n  =  0;  n  <  count [cliqueNumber  +  1];  n++)  { 

System. out .print ( (n  +  1)  +  ")  "); 

for  (int  m  =  0;  m  <  cliqueNumber  +  1;  m++)  { 

System. out. print ( (clique [cliqueNumber  +  1] [n] [m] 

max imumCli quel [m]  =  clique [cliqueNumber  + 


System. out. print In () ; 

} 

System. out. print In () ; 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
k  k  k 


// 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


k  kkkk  k 


+  1) 
k)  )  ; 


n++) 


{ 

1)  - 


I? 


k] 


for  (int  k  =  0;  k  <  cliqueNumber;  k++)  { 

System. out . println ( "There  are  "  +  count [ (cliqueNumber 

+  "  Cliques  K"  +  ((cliqueNumber  +  1)  - 


System. out. println ("The  Cliques  are:  "); 

//  for  (int  k  =  0;  k  <vectorsNumber ;  k++)  { 

for  (int  n  =  0;  n  <  count [ (cliqueNumber  +  1)  -  k]  ; 

System. out. print ( (n  +  1)  +  ")  "); 

for  (int  m  =  0;  m  <  (cliqueNumber  +  1)  -  k;  m++) 


k]  [n]  [m]  +  1) 


System. out .print ( (clique [ (cliqueNumber  + 

+  "  " )  ; 

//  System. out . print ( (clique [ k] [n] [m] +1 )  + 


} 

System. out. print ln() ; 

} 

System. out .println ( ) ; 

} 


} 


} 


49 


THE  OUTPUT 


Maximum  Clique  Problem  of  G(37,85) 
Find  the  degree  of  each  vertex 


degree 

(vl  )  = 

=  2 

degree 

(v2  )  = 

=  4 

degree 

(v3  )  = 

=  4 

degree 

(v4  )  = 

=  3 

degree 

(v5  )  = 

=  4 

degree 

(v6  )  = 

=  2 

degree 

(v7  )  = 

=  2 

degree 

(v8  )  = 

=  6 

degree 

(v9  )  = 

=  3 

degree 

(vlO  ) 

=  4 

degree 

(vll  ) 

=  3 

degree 

(vl2  ) 

=  6 

degree 

(vl3  ) 

=  3 

degree 

(vl4  ) 

=  3 

degree 

(vl5  ) 

=  3 

degree 

(vl 6  ) 

=  1 

degree 

(vl7  ) 

=  14 

degree 

(vl  8  ) 

=  5 

degree 

(vl 9  ) 

=  3 

degree 

(v20  ) 

=  10 

degree 

(v21  ) 

=  7 

degree 

(v22  ) 

=  10 

degree 

(v23  ) 

=  15 

degree 

(v24  ) 

=  3 

degree 

(v25  ) 

=  4 

degree 

(v2  6  ) 

=  4 

degree 

(v27  ) 

=  2 

degree 

(v2  8  ) 

=  8 

degree 

(v29  ) 

=  5 

degree 

(v30  ) 

=  10 

degree 

(v31  ) 

=  3 

degree 

(v32  ) 

=  4 

degree 

(v33  ) 

=  1 

degree 

(v34  ) 

=  1 

degree 

(v35  ) 

=  4 

degree 

(v36  ) 

=  2 

degree 

(v37  ) 

=  2 

Potential  Max  Clique  Degree:  16. 
There  is  not  a  K16  Max  Clique. 
Check  for  K15  Max  Clique. 

There  is  not  a  K15  Max  Clique. 
Check  for  K14  Max  Clique. 

There  is  not  a  K14  Max  Clique. 
Check  for  K13  Max  Clique. 

There  is  not  a  K13  Max  Clique. 
Check  for  K12  Max  Clique. 
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There  is  not  a  K12  Max  Clique. 
Check  for  Kll  Max  Clique. 

There  is  not  a  Kll  Max  Clique. 
Check  for  K10  Max  Clique. 

There  is  not  a  K10  Max  Clique. 
Check  for  K9  Max  Clique. 

There  is  not  a  K9  Max  Clique. 
Check  for  K8  Max  Clique. 

There  is  not  a  K8  Max  Clique. 
Check  for  K7  Max  Clique. 

Maybe  there  is  a  K7  Max  Clique. 

There  is  a  K2  Clique. 

The  Clique  is: 

1  2 

There  is  a  K3  Clique. 

The  Clique  is: 

12  5 


The  Maximum  Clique  is  K6 
There  is  1  Maximum  Clique 
The  Maximum  Clique  is: 


1) 

17 

18 

21 

22 

23 

28 

There 

is 

1  ( 

Clique 

K6 

The 

Cliques 

are : 

1) 

17 

18 

21 

22 

23 

28 

There 

are  9 

Cliques  K5 

The 

Cliques 

are : 

1) 

17 

18 

21 

22 

23 

2) 

17 

18 

21 

22 

28 

3) 

17 

18 

21 

23 

28 

4) 

17 

18 

22 

23 

28 

5) 

17 

21 

22 

23 

28 

6) 

17 

21 

23 

25 

28 

7) 

17 

22 

23 

28 

32 

8) 

17 

22 

23 

29 

30 

9) 

18 

21 

22 

23 

28 

There  are  31  Cliques  K4 
The  Cliques  are: 

1)  8  12  13  20 

2)  15  17  21  23 

3)  17  18  21  22 

4)  17  18  21  23 

5)  17  18  21  28 
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6) 

17 

18 

22 

23 

7) 

17 

18 

22 

28 

8) 

17 

18 

23 

28 

9) 

17 

21 

22 

23 

10) 

17 

21 

22 

28 

11) 

17 

21 

23 

25 

12) 

17 

21 

23 

28 

13) 

17 

21 

25 

28 

14) 

17 

22 

23 

28 

15) 

17 

22 

23 

29 

16) 

17 

22 

23 

30 

17) 

17 

22 

23 

32 

18) 

17 

22 

28 

32 

19) 

17 

22 

29 

30 

20) 

17 

23 

25 

28 

21) 

17 

23 

28 

32 

22) 

17 

23 

29 

30 

23) 

18 

21 

22 

23 

24) 

18 

21 

22 

28 

25) 

18 

21 

23 

28 

26) 

18 

22 

23 

28 

27) 

20 

24 

26 

31 

28) 

21 

22 

23 

28 

29) 

21 

23 

25 

28 

30) 

22 

23 

28 

32 

31) 

22 

23 

29 

30 

There 

are 

66 

Cliques  K3 

The 

Cliques 

are : 

1) 

1  2 

5 

2) 

2  3 

4 

3) 

3  4 

8 

4) 

8  12  13 

5) 

8  12  14 

6) 

8  12  20 

7) 

8  13  20 

8) 

9  10  17 

9) 

10 

17 

23 

10) 

11 

17 

23 

11) 

12 

13 

20 

12) 

15 

17 

21 

13) 

15 

17 

23 

14) 

15 

21 

23 

15) 

17 

18 

21 

16) 

17 

18 

22 

17) 

17 

18 

23 

18) 

17 

18 

28 

19) 

17 

19 

22 

20) 

17 

21 

22 

21) 

17 

21 

23 

22) 

17 

21 

25 

23) 

17 

21 

28 

24) 

17 

22 

23 

25) 

17 

22 

28 

26) 

17 

22 

29 

52 


27) 

17 

22 

30 

28) 

17 

22 

32 

29) 

17 

23 

25 
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APPENDIX  E.  THE  CODE  OF  THE  PRUNING  ALGORITHM 


/*  *  File:  Maximum  Clique 

* 

*  Created  on  March  03,  2008 

* 

*  Modified  on  March  012,  2008 

* 

* 

*  ***********THIS  WORKS  FOR  the  GRAPH  "911",  G(37,85)  where  there 
is  a  K6  (17,18,21,22,23,28) 

* 

* 

*  Let  G= (V (G) , E (G) )  be  a  graph  with  n=|V(G)|  vertices  and 
m= | E (G) I  edges . 

*  Let  v  i  belongs  to  V(G) .  Then  N [V  i]=set  of  neighbors  of  v  i, 
together  with  v  i  itself,  called  the  closed  neighborhood.  The  subgraph 
induced  by*  N[V_i]  is  denoted  by  <N[V  i]>. 

*/ 

public  static  void  main ( String [ ]  args)  { 

//  initialize  instance  variables 
int  vectorsNumber  =  37; 

int [ ]  degree  =  new  int [vectorsNumber ] ; 

int [ ]  booleanDegree  =  new  int [vectorsNumber ] ; 

int  maxCliqueDegree  =  0; 

int  tempMaxCliqueDegree  =  0;//  This  is  equal  to  (Potential 
Maximum  Clique  #  -1),  (PMC#-1) .  Its  value  decreases  upon  the  condition 

which  are  met,  until  Maximum  Clique  is  found.  Then  Clique#= (updated  by 
the  algorithm) PMC# 

int  initialTempMaxCliqueDegree  =  0;//  This  is  equal  to  the 
initial  value  of  the  (Potential  Maximum  Clique  #  -1),  (PMC#-1)  .  Its 

value  remains  constant. 

int  newTempMaxCliqueDegree  =  0; 
boolean  maxClique  =  false; 
boolean  initialMaxClique  =  false; 

int[][]  table  =  new  int [vectorsNumber] [vectorsNumber] ; // 
This  is  the  nxn  adjacency  matrix  of  the  given  graph  G  (with  n  vertices) . 
int [ ] []  maximalClique  =  new  int [1000] [vectorsNumber] ; 
int  countMaximalCliques=0 ; 
int  row,  column; 

//********************************************************************** 

'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k-k-k-k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k'k'k'k'k'k'k 

//  I.  PUT  AN  UPPER  BOUND 

//  1.  Set  the  nxn  adjacency  matrix  of  the  given  graph  G 

(with  n  vertices) . 

table  [0]  [0]  =  0; 
table  [0]  [1]  =  1; 
table  [0]  [2]  =  0; 
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table [36] [34]  =  1; 
table  [36]  [35]  =  0; 
table [36] [36]  =  0; 

II  **************** 


System. out . println ( "Maximum  Clique  Problem  of  G(37,85)"); 
System. out. print In () ; 

//  2.  Find  the  degree  of  each  vertex  by  adding  the  entries 
of  each  row  of  the  adjacency  matrix. 

System. out . println ( "Find  the  degree  of  each  vertex  "); 
for  (row  =  0;  row  <  vectorsNumber ;  row++)  { 

int  sum  =  0;//  This  is  the  degree  of  each  vertex 
for  (column  =  0;  column  <  vectorsNumber;  column++) 
sum  =  sum  +  table [row] [column] ; 
degree [row]  =  sum; 

} 

//  print  the  degree  of  each  vertex 

for  (row  =  0;  row  <  vectorsNumber;  row++)  { 

System. out. println ("degree (v"  +  (row  +  1)  +  "  )=  "  + 

degree [row] ) ; 

} 

System. out. print ln() ; 

//  3.  Find  the  Potential  Max  Clique#,  say  PMQ#, 

//  (Potential  Max  Clique#  =  max  Vertex  Degree+1) . 
for  (row  =  0;  row  <  vectorsNumber;  row++) 

if  (degree [row]  >  tempMaxCliqueDegree)  { 
tempMaxCliqueDegree  =  degree [row] ; 

} 

System. out . println (" Potential  Max  Clique  Degree:  " 

+  (tempMaxCliqueDegree  +  1)  + 


initialTempMaxCliqueDegree  =  tempMaxCliqueDegree; 

//  4 .  Let  v  k  the  vertices  with  degree  greater  than  or  equal 
to  (PMC#-l),for  some  k,  l<=k<=n.  Find  the  #  of  vertices  v  k  with  degree 
greater  than  or  equal  to  (PMC#-1),  say  "inducedSum" . 

do  { 

int  inducedSum  =  0;//  the  #  of  vertices  with  degree 
greater  than  or  equal  to  (PMC#-1) 

for  (row  =  0;  row  <  vectorsNumber;  row++)  { 

if  (degree [row]  >=  tempMaxCliqueDegree)  { 
booleanDegree [row]  =  1; 

inducedSum  =  inducedSum  + 

booleanDegree [row] ; 

}  else 

booleanDegree [row]  =  0; 

} 

//  4a.  If  "inducedSum"  is  greater  than  or  equal  to 

PMC#,  then  maybe  there  is  a  Max  Clique  with  Clique#  =  PMC#. 

if  (inducedSum  >=  tempMaxCliqueDegree  +1)  { 

System. out .println ( "Maybe  there  is  a  K" 
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+  (tempMaxCliqueDegree  +1)  +  "  Max  Clique."); 
System. out .println ( ) ; 
initialMaxClique  =  true; 

//  4b.  Else  if  "inducedSum"  is  NOT  greater  than 
or  equal  to  PMC#,  there  is  not  a  Maximum  Clique  with 
Clique#  =  PMQ# .  Set  PMQ#=PMQ#-1  and  continue  by  going 
back  to  step  4 . 

}  else  { 

System. out . println ( "There  is  not  a  K" 

+  (tempMaxCliqueDegree  +1)  +  "  Max  Clique."); 
System. out .println ( "Check  for  K" 

+  tempMaxCliqueDegree)!  "  Max  Clique."); 

System. out. print ln() ; 
tempMaxCliqueDegree--  ; 

} 

}  while  (tempMaxCliqueDegree  >  0  &&  initialMaxClique  == 

false) ; 

// 

************************************************************************ 

'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

//  II.  Search  for  the  Max  Clique,  based  on  the  Upper  Bound 

newTempMaxCliqueDegree  =  tempMaxCliqueDegree;//  This  is 
equal  to  (Potential  Maximum  Clique  #  -1),  (PMC#-1) . 

//  Its  value  decreases  at  the  end  of  the  algorithm  if  all 
the  intermediate  conditions  are  NOT  met,  until  Maximum  Clique  is  found. 
//  Then  Clique#= (updated  by  the  algorithm) PMC# 

int [  ]  adj acentVector  =  new  int [ initialTempMaxCliqueDegree  + 
1];//  The  vector  where  the  vertices  of  the  closed  neighborhood  of  each 
v  k  are  stored. 

int [ ]  booleanDegreeOf InducedTable  =  new 

int [ initialTempMaxCliqueDegree  +  1]; 

int [ ]  degreeOf InducedTable  =  new 

int [ initialTempMaxCliqueDegree  +  1]; 

//  1.  For  every  v  k  with  degv  k=(PMC#-l),  investigate  if 

<N[Vk]>  is  a  Max  Clique  with  Clique#  =  PMQ#: 

do  { 

System. out. println ("Check  for  K"  + 

(newTempMaxCliqueDegree  +  1) 

+  "  Max  Clique."); 

System. out .println ( ) ; 
int  countForMaxClique  =  0; 
int  vectorOfMaxClique  =  0; 

//  a.  Build  <N[V  k]>: 

do  { 

int  count  =  0;//  This  is  the  cardinality  of 

v  k '  s 

//  neighborhood  set. 

vectorOfMaxClique  =  countForMaxClique; 
adj acentVector [ 0 ]  =  countForMaxClique; 
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System. out .println ( ) ; 

System. out .println ( "Check  if  the  vertex  V" 

+  (countForMaxClique  +1)  +  "  belongs  to  a  K" 

+  (newTempMaxCliqueDegree  +1)  +  "  Max  Clique."); 

adj acentVector [count]  =  countForMaxClique; 

//  i.  First,  find  the  adjacent  vertices  to  v  k. 
//  ii.  Then,  store  v  k  and  its  adjacent  vertices 
in  a  one  dimensional  array,  say  [V_k]  (adj acentVector [] ) . 

for  (column  =  0;  column  <  vectorsNumber ; 

column++)  { 

if  (table [vectorOfMaxClique] [column]  ==  1) 

{ 

count++; 

adj acentVector [count]  =  column; 

} 

} 


//  b.  Check  if  the  degrees  of  the  vertices  of 
<N [V  k] >  are  greater  than  or  equal  to  the  (PMQ#-1) : 

//  i.  First,  find  the  degree  in  <N [V  k] >  of  each 

vertex  of  <N [V  k]> 

for  (row  =  0;  row  <  count  +  1;  row++)  { 

int  suml  =  0;//  This  is  the  degree  in 
<N [V  k] >  of  each  vertex  of  <N [V  k] > 

for  (column  =  0;  column  <  count  +  1; 


column++) 


suml  =  suml  + 

table [adj  acentVector [row] ]  [adj  acentVector [column] ] ; 

degreeOf InducedTable [row] 


} 


suml  ; 


//  ii.  Then,  find  the  #  of  vertices  with  degree 
greater  than  or  equal  to  the  (PMQ#-1) ,  say  "inducedSuml " . 

int  inducedSuml  =  0;//This  is  the  #  of  vertices, 
of  <N[V^k]>  ,  with  degree  greater  than  or  equal  to  the  (PMQ#-1) . 

for  (row  =  0;  row  <  count  +  1;  row++)  { 

if  (degreeOf InducedTable [row]  >= 

newTempMaxCliqueDegree)  { 

booleanDegreeOf InducedTable [row]  = 

1; 

inducedSuml ++; 

}  else 

booleanDegreeOf InducedTable [row]  = 

0; 

/ /  inducedSuml  =  inducedSuml  + 

/ /  booleanDegreeOf InducedTable [row] ; 

} 

//  c.  If  the  "inducedSuml"  is  greater  than  or 
equal  to  the  PMQ#  then  the  Clique#  =  PMQ# . 

if  (inducedSuml  >=  newTempMaxCliqueDegree  +1)  { 
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System. out . println ( "YES !! !  V"  + 

(countForMaxClique  +  1)  +  "  belongs  to  at  least  one  K" 
+ (newTempMaxCliqueDegree  +1)  +  "  Max  Clique."); 

System. out. print In ( ) ; 

int  rowUpTriangular  =  1; 
int  columnUpTriangular  =  1; 

//  Find  the  Maximum  Cliques  that  v  k 
belongs  to.  Consider  the  one  dimensional  array  [V  k]  , 

( "adj acentVector [ ] " ) ,  which  contains  the  elements  of  the  closed 

neighborhood  of  v  k,  N[V  k] . 

/  /  *  *  *  *  *  *  * 

//  Also,  consider  each  element  of 
v  k' s  neighborhood,  say  v  m,  for  some  m,  l<=m<= (cardinality  of  v  k's 

//  neighborhood  set) . 

for  (row  =  rowUpTriangular;  row  < 

count  +  1;  row++)  { 

int  top  =  -1;//  counts  the  top  for 
the  vector  with  the  vertices  adjacent  to  the  "second  vertex" 

int  topMaxClique  =  2;// 

counts  the  top  for  the  vector  with  the  Max  Clique  vertices 

int[]  stack  =  new 
int [ initialTempMaxCliqueDegree  +  1];//  This  is  the  vector  where  the 

vertices  of  the  neighborhood  of  each  element  of  N [V  k]  are  stored. 

int [ ]  stackMaxClique  = 

new  int  [  initialTempMaxCliqueDegree  +  1];//  This  is  the  vector  where  the 
vertices  consisting  a  Maximum  Clique  are  stored. 

System. out . print ln() ; 
System. out . print In ( "the 

element  v  m  now  is:  "  +  (adj acentVector [row]  +  1) ) ; 

System. out. print ln() ; 


//  1)  For  each  element 

of  v  k's  neighborhood, say  v  m,  for  some  m,  l<=m<= (cardinality  of  v  k's 
neighborhood  set),  find  its  neighborhood  N (V  m)  in  <N [V  k]>  and  store  it 
in  a  one  dimensional  array,  say  "stack". 

for  (column  = 

columnUpTriangular;  column  <  count  +  1;  column++)  { 

if 

(table [adj acentVector [row] ] [adj acentVector [column] ]  ==  1)  { 


top++; 

stack 


adj acentVector [column]  ; 

System. out . println ( "stack [ "  +  top  +  "]=  "  +  stack 

} 


System. out . println ( "top=  "  + 


} 


System. out. println() 


//  2)  Create  another  one 

dimensional  array,  say  "stackMaxClique".  This  is  the  vector  where  the 
vertices,  consisting  a  Maximum  Clique,  are  stored.  Store  v  k  as  first 
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element,  v  m  as  second  element  (starting  with  m=l)  and  the  first  element 
of  the  "stack",  as  third  element  (these  3  vertices  consist  a  K3) . 

stackMaxClique [ 0 ]  = 


countForMaxClique; 
adj acentVector [row] ; 


stackMaxClique [1] 
stackMaxClique [2 ] 


stack [ 0 ] ; 


//  3)  Set  "elementA"=the  first  element  of  the  "stack",  and 

"elementB"=the  second  one.  Check  if  "elementA"  and  "elementB"  are 
adj  acent : 

int  countTop  =  0; 
int  stackElementA  =  0; 
int  stackElementB  =  1; 
boolean  stop  =  false; 


if  (top  <  1)  { 


System. out . println ( "For  v  k=  "  +  (countForMaxClique  +1)  +  "  and 

'  +  (adj acentVector [row]  +1)  +  "  there  is  NOT  a  Max  Clique."); 

} 


if  (top  >=  1)  { 

do  { 

if 

(table [ stack [ stackElementA] ][ stack [ stackElementB] ]  ==  1)  { 


countlnTheStackMaxClique  =  0; 


int 


//  If  yes,  check  if  the  "elementB"  is  adjacent 
"stackMaxClique"  (start  checking  by  the 
"stackMaxClique") . 

(int  i  =  0;  i  <=  topMaxClique;  i++)  { 


//  a) 

with  each  element  of 
third  element  of 

for 


if  (table [stackMaxClique [i] ] [stack [stackElementB] ]  ==  1)  { 


countlnTheStackMaxClique+t; 


} 


} 


System. out . println ( "topMaxClique=  "  +  topMaxClique); 

System. out. print In ( "countInTheStackMaxClique=  " 

+  countlnTheStackMaxClique) ; 

//  i.  If  yes,  then  add  "elementB"  to  "stackMaxClique",  and  go  back  to 
step  4),  where  "elementA"="elementB"  and  "elementB"=  the  next  element  of 
"stack". 
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if 

(countlnTheStackMaxClique  ==  topMaxClique  +1)  { 

topMaxClique++; 

stackMaxClique [topMaxClique]  =  stack [ stackElementB] ; 

stackElementA  =  stackElementB; 

stackElementB++; 

System. out . println (" stackMaxClique [ "  +  topMaxClique  +  "]=  " 

+  stackMaxClique [topMaxClique]  )  ; 


countTop++; 


} 


//  ii.  If  no  then  go  back  to  step  4),  where  "elementA"  remains  the  same 
and  "elementB"=the  next  element  of  "stack". 


stackElementB++; 

//  This  is  the  condition  to  terminate  the  do-while  loop. 

((stackElementA  ==  top  &&  stackElementB  ==  (top  +  1)) 

I |  (stackElementA  ==  (top  -  1)  &&  stackElementB  ==  (top  +  1))) 


else  { 

} 

if 


stop  =  true; 


} 


//  b)  If  no  go  back  to  step  4),  where  "elementA"  remains  the  same  and 
"elementB"=the  next  element  of  "stack". 

else  { 


stackElementB++; 


} 


//  c)  If 

check  all  elements  of  "stack",  as  "elementB"  then  "elementA"=the  next 
element  of  "stack"  and  "elementB"=the  element  of  "stack"  which  is  after 
the  new  "elementA". 

if 

(stackElementB  ==  top  +1)  { 

stackElementA++; 

stackElementB  =  stackElementA  +  1; 

} 

System. out . println ( "stop=  "  +  stop) ; 

//  This  is  the  condition  to  terminate  the  do-while  loop. 


if 

((stackElementA  ==  top  &&  stackElementB  ==  (top  +  1)) 

I |  (stackElementA  ==  (top  -  1)  &&  stackElementB  ==  (top  +1)))  { 


true  ; 


} 


stop  = 


System. out . println ( "stackElementA=  "  +  stackElementA); 
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System. out. println ( "stackElementB: 


+  stackElementB) ; 


false) ; 


}  while  (stop  == 

if (topMaxClique>maxCliqueDegree) { 
maxCliqueDegree=topMaxClique+l ; 

} 

for  (int  j  =  0;  j  <=  topMaxClique;  j++)  { 

maximalClique [countMaximalCliques] [ j ] = 
stackMaxClique [ j ]  +  1; 

} 

countMaximalCliques++; 


<= 


of 


System. out . print ( "The  Max  Clique  is:  "); 
topMaxClique;  i++)  { 

System. out . print ( "V"  +  ( stackMaxClique [ i ] 

} 


for  (int  i  = 

+  1 )  +  "  " )  ; 

} 


0;  i 


} 

// 


N [V  k] ,  v  (m+1),  going  back  to  step 

} 


System. out. print In () ; 

} 


d)  Continue  with  the  next  element 
1)  • 


//********************************************************************** 

******************************* 

//  v.  Else  if  the  "inducedSuml"  is  equal 
to  the  PMQ#  then  v  k  belongs  to  one  Max  Clique  with  Clique#  =  PMQ# . 

else  if  (inducedSuml  == 

newTempMaxCliqueDegree  +1)  { 

System. out. println ("YES !! !  There  is 
a  K"  +  (newTempMaxCliqueDegree  +1)  +  "  Max  Clique."); 

maxClique  =  true; 

maxCliqueDegree  = 

newTempMaxCliqueDegree  +  1; 

System. out . print ( "The  Max  Clique  is: 

V"  +  (countForMaxClique  +1)  +  "  "); 


newTempMaxCliqueDegree)  { 

(adj acentVector [i]  +  1)  +  "  "); 

} 

} 


for  (int  i  =  1;  i  <=  count;  i++) 

if  (degreeOf InducedTable [ i ]  >= 

System. out . print ( "V"  + 


} 

System. out. print ln() ; 
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//  d.  If  the  #  of  vertices  of  <N[V^k]>  are  not 
greater  than  or  equal  to  the  PMQ#  then  continue  searching  for  a  Clique# 
=  PMQ#  by  going  back  to  step  II. 1.  with  the  next  vertex  v  (k+1) 

else 

System. out .println ( "The  vertex  V"  + 
(countForMaxClique  +  1)  "  does  not  belong  to  a  K" 

(newTempMaxCliqueDegree  +1)  +  "  Max  Clique.  "); 

count ForMaxClique++; 

}  while  (countForMaxClique  <  vectorsNumber ) ; 

System. out . print In () ; 

//  II.  2.  If  there  is  not  a  Clique#  =  PMQ#  then  set 
PMQ#  =  PMQ#-1  and  continue  searching  for  a  Maximum  Clique  with  Clique#  = 
(PMC#-1)  by  going  back  to  step  1.4. 

newTempMaxCliqueDegree-- ; 

}  while  (newTempMaxCliqueDegree  >  0  &&  maxClique  ==  false)  ; 

/ /  print  the  outputs 

System. out . println ( "The  Maximal  Cliques  are:"); 
for  (int  i  =  0;  i  <  countMaximalCliques ;  i++)  { 

System. out .println ( )  ; 

for  (int  j  =  0;  j  <  maxCliqueDegree;  j++)  { 

System. out .print (maximalClique [i] [j]+"  "); 

} 

} 

System. out .println ( ) ; 

System. out . println ( "The  Maximum  Clique  is  K"  + 
maxCliqueDegree) ; 

} 

} 
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THE  OUTPUT 


Maximum  Clique  Problem  of  G(37,85) 
Find  the  degree  of  each  vertex 


degree 

(vl  )  = 

=  2 

degree 

(v2  )  = 

=  4 

degree 

(v3  )  = 

=  4 

degree 

(v4  )  = 

=  3 

degree 

(v5  )  = 

=  4 

degree 

(v6  )  = 

=  2 

degree 

(v7  )  = 

=  2 

degree 

(v8  )  = 

=  6 

degree 

(v9  )  = 

=  3 

degree 

(vlO  ) 

=  4 

degree 

(vll  ) 

=  3 

degree 

(vl2  ) 

=  6 

degree 

(vl3  ) 

=  3 

degree 

(vl4  ) 

=  3 

degree 

(vl5  ) 

=  3 

degree 

(vl 6  ) 

=  1 

degree 

(vl7  ) 

=  14 

degree 

(vl  8  ) 

=  5 

degree 

(vl 9  ) 

=  3 

degree 

(v20  ) 

=  10 

degree 

(v21  ) 

=  7 

degree 

(v22  ) 

=  10 

degree 

(v23  ) 

=  15 

degree 

(v24  ) 

=  3 

degree 

(v25  ) 

=  4 

degree 

(v2  6  ) 

=  4 

degree 

(v27  ) 

=  2 

degree 

(v2  8  ) 

=  8 

degree 

(v29  ) 

=  5 

degree 

(v30  ) 

=  10 

degree 

(v31  ) 

=  3 

degree 

(v32  ) 

=  4 

degree 

(v33  ) 

=  1 

degree 

(v34  ) 

=  1 

degree 

(v35  ) 

=  4 

degree 

(v36  ) 

=  2 

degree 

(v37  ) 

=  2 

Potential  Max  Clique  Degree:  16. 
There  is  not  a  K16  Max  Clique. 
Check  for  K15  Max  Clique. 

There  is  not  a  K15  Max  Clique. 
Check  for  K14  Max  Clique. 

There  is  not  a  K14  Max  Clique. 
Check  for  K13  Max  Clique. 

There  is  not  a  K13  Max  Clique. 
Check  for  K12  Max  Clique. 
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There  is  not  a  K12  Max  Clique. 

Check  for  Kll  Max  Clique. 

There  is  not  a  Kll  Max  Clique. 

Check  for  K10  Max  Clique. 

There  is  not  a  K10  Max  Clique. 

Check  for  K9  Max  Clique. 

There  is  not  a  K9  Max  Clique. 

Check  for  K8  Max  Clique. 

There  is  not  a  K8  Max  Clique. 

Check  for  K7  Max  Clique. 

Maybe  there  is  a  K7  Max  Clique. 

Check  for  K7  Max  Clique. 

Check  if  the  vertex  VI  belongs  to  a  K7  Max  Clique. 
The  vertex  VI  does  not  belong  to  a  K7  Max  Clique. 

Check  if  the  vertex  V2  belongs  to  a  K7  Max  Clique. 


The  vertex  V37  does  not  belong  to  a  K7  Max  Clique. 
Check  for  K6  Max  Clique. 

Check  if  the  vertex  VI  belongs  to  a  K6  Max  Clique. 
The  vertex  VI  does  not  belong  to  a  K6  Max  Clique. 

Check  if  the  vertex  V2  belongs  to  a  K6  Max  Clique. 


Check  if  the  vertex  V17  belongs  to  a  K6  Max  Clique. 
YES! ! !  There  is  a  K6  Max  Clique. 

The  Max  Clique  is:  V17  V18  V21  V22  V23  V28 


The  Maximum  Clique  is  K6 
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